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The Theory of Plasticity in Relation to its 
Engineering Applications 
f : By Hueu Forp, London!) 


Nature of Plastic Deformation 


There are two main approaches to the study of plastic deformation: these 
_may be termed the ‘microscopic’ and the ‘macroscopic’. 

The ‘microscopic’ line of study is that which is variously called ‘metal 

physics’, ‘physics of single crystals’, etc. and which relates the behaviour of 
_ crystalline aggregates to fundamental ideas of the deformation of single crystals. 

It is possible to formulate some general notions of the theoretical strength and 

cohesion of materials and of the nature of shear arising from dislocations in 
regular crystal structure [1], [2], [3]*). Lately Hirt and BisHop [4] have been 

preparing the way for building up a mathematical theory of plasticity in poly- 

crystals from single crystal behaviour, but so far only general theorems are 
_ being studied. 

The ‘macroscopic’ approach is that in which the material is assumed to be 
an agglomerate without a crystalline structure. The mathematical theory of 
plasticity deals with idealized materials of this kind, and examines analytically 
their behaviour when subjected to complex stress conditions. Its purpose is to 
develop.mathematical techniques for calculating non-uniform distributions of 
stress and strain in bodies permanently distorted in any way. Leading out of 

the strict theory is the ‘technological’ study: use is made on the one hand of 
' the relationships of the mathematical theory and on the other of experimentally 
determined laws and criteria to establish technically important results. 

A full treatment cannot be given in one paper: all that can be done is to 
discuss the basic laws of the mathematical theory and review some of the 
techniques used for solving problems. A few examples of the application of 
these laws and techniques to the study of industrial processes will then be given. 

One of the difficulties in interpreting the mathematical theory of plasticity 

_is that the behaviour of technical materials when under load in the plastic range 
depends upon their structure and the problems presented by a complete ana- 
lysis of anisotropic, non-homogeneous materials are too complex to attempt a 
solution by our present methods. It is in fact remarkable that metallic materials 
behave so nearly to that predicted by theory from consideration of idealized 


1) Department of Applied Mechanics, Imperial College, University of London. 
2) Numbers in brackets refer to References, page 34. 
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materials, but we must always remember that the anisotropic, crystalline nature 
of metals must intrude sooner or later and modify the results predicted from — 
simple theory. ; : 
The mathematical theory must, therefore, make certain assumptions and 
generally these are:— | 
(1) Isotropy. Theory for anisotropic material is not yet very far advanced [5]. | 
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(2) Sharp yield with no hysteresis on loading and unloading. No change in 
elastic moduli with the degree of deformation. 

(3) No Bauschinger effect —one basic yield stress relation for tension and com- 
pression. 

(4) No temperature effects. Yield stress independent of time rate of strain, no 
creep or recovery. 

(5) Idealized stress-strain relations (Figure 1). 

In large deformations it is possible to obtain complete solutions by a step- 
wise graphical solution, taking small strain increments over which the yield 
stress is considered constant (Figure 2). 

Actual stress-strain curves can be approximated to types (c) and (d) with 
reasonable success and although the sharp yield does not usually occur, suitable 
experimental techniques can be used, as will be shown later. 

The theory of perfectly plastic solids with no elastic region can be satisfac- 
torily applied to work-hardening materials. This theory has wide applications, 
particularly for large deformations where the elastic strains are an insignificant 
part of the whole. When we determine that, even for the highest strength steel, 


the elastic strain does not exceed 0-5 per cent, this theory applies very generally 
in practical problems. 
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There are two principal types of elastic-plastic problem: That in which all 
the material reaches yield at the same instant, but yielding depends upon the 
_ elastic strain pattern, the plastic strains being of the same magnitude, and that 
in which there is an elastic-plastic interface and the elastic region presents a 
restraint to unrestricted flow. The equations of flow involve the generalized 
Hooke’s law equations, and the elastic and plastic strains are of the same 


Stress 


Strain 


Figure 2 


order. Only a few cases have been solved, a high degree of symmetry usually 
being necessary, and the shape of the elastic-plastic interface must be pre- 
dictable [5]. 


The General State of Stress 


The stress at a point can be defined in terms of six stress components. 
Three normal stresses 


Paws Puy Pus 


and three shear stresses 


Pay, Puz Pew 


according to any arbitrary Cartesian co-ordinates [6]. 
Alternatively, a state of stress is uniquely defined by its three principal 
stresses: 


Pr Px» Ps 


and the directions of these stresses in space. 

A given state of stress at a point does not depend upon the directions of the 
arbitrary co-ordinates chosen in the first place, and, therefore, there can only 
be one set of principal stresses and axes. 

In terms of the general stress components poe, Pyy «++» Pza the principal 
stresses are the three roots of the equation 


P+hrt+hp+ss=0 (1) 
where 


h= — (Peo + Puy + Pes) » (la) 
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Js ss (Pew Pyy + Puy Paz + Pus Pea — piv =~ Piz es Dix) , (1b) 
J3 we (Pow Puy Paz Sip 2 Pay Pye Pie a Pow Pez sw Puy Piz >= Paz Ms (1c) 


For a given state of stress the quantities f,, J,, and J, must be constant, i.e. 
they are invariants to orthogonal transformations of co-ordinates. 
Notation: pz, = stress in direction x on a plane x (normal), 
bey = stress in direction x on a plane y (shear). 
Principal stresses p, > p. > p; algebraically. 
Tension is positive, compression negative, unless otherwise indicated. 


Tensor Notation 


In the generalized statement of the mathematical theory of plasticity, tensor 
notation allows a clear and concise presentation, but, since the greater part of 
the paper deals with particular problems and processes, it will not be used here. 
Reference should be made to Hitt [5]. 


Yield 


Suppose the intensity of a state of stress at a point in an idealized material 
gradually increases, passing through a series of stress states, when does the 
first permanent change of shape occur? 

For isotropic materials, over the range of stresses normally encountered in 
technology, yielding is found to be independent, over the practically significant 
range, of hydrostatic stress (i.e. the mean stress) being dependent only on 
stress differences [7]. 

Denote the hydrostatic stress by 


f= Ramee Gaesh Cael yest sas, (2) 
The stress deviations are then 
. Se god) ea 
Ds —— Dux pb ee (Deve oe Bea os (3) 
Note: 
Pay + Pay : 
In terms of the stress deviations (3), equation (1) becomes 
Po + ip +h=o, (4) 
where 
Je = (boa buy + buy Pie + bie Pow — P2y — P22 — Phe)» (4b) 


Js =F (Pave Pav Pas +2 Pey Pus Puw mr Poe Pin Pr Pay Pi. = Dis Pay) c (4c) 
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Now yielding can only depend on the magnitude of the state of stress and 
must be independent of the orientation of the co-ordinate axes, i.e. it can only 


_ bea function of Jz and J,. 


Also we have assumed that there is no Bauschinger effect so yield must not 
depend upon change of sign of the stresses. 

Jz is unaltered by change of sign, but Jz is, and yielding would, therefore, 
need to be an even function of Jj. This would soon lead to complicated func- 
tions for yielding. While it must be remembered that yielding could be any 
suitable function of Jz and Jz the simplest is clearly to take 


Jz = const. 


This is the von Mises yield criterion; CLERK-MAXWELL in 1865 first pro- 
posed it on mathematical grounds, but other interpretations have been put 
forward [8], [9], [10]. 

By carrying out a simple tensile or compression test, the yield stress can be 
determined. Alternatively, torsion of a thin tube allows the shear yield stress 
in torsion to be found. 

Denoting the shear yield stress in the tension or compression test by k 
(i.e. the yield stress measured in the test would be 2 k) the von Mises criterion is 


meen SS bx — Ba)? + ba — Pa) + (be — Bi)? = 8 (5) 
or, for general co-ordinates as: 

(Paw — Pyy)* + (buy — Bez)? + (Pez — Puw)* + © Pay + 6 Diz + 6 D7, = 8 k*. (Sa) 
Based on the torsion shear stress k*, we obtain: 


(1 — bo)? + (be — bs)” + (b3 — f1)? = 6 R*?. (5b) 


The other practically important criterion is the Tvesca criterion, sometimes 
called GuEstT’s law because GuEsT (1900) first appreciated its practical value 
and carried out accurate tests with combined stress conditions. It is practically 
simpler and takes od: 


This criterion ignores the effect of the intermediate stress f,. The greatest 
difference between the two criteria is about 15-5 per cent. For a great number 
of problems, particularly for plane strain, the two criteria give the same result 
except for a constant. The constant can be adjusted so that results never differ 
by more than 7-75 per cent. 

For rigid-plastic or elastic-plastic materials we take J, = const for all stress 
states irrespective of total strain. For work-hardening materials, J, = const 
only for an increment of strain. The constant then becomes a parameter, deter- 
minable from the yield stress-strain curve, in a manner to be explained later. 
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Limiting Stress Fields 


PRAGER and HopcE [11] and Ht [12] have discussed the graphical repre- 
sentation of the limiting conditions of yield by the Tresca and von Mises criteria. 
These representations are useful in considerations of testing plasticity laws. 


The Law of Plastic Flow 


When yield has occurred, what is the resulting strain distribution ? 

If the stress is uniform, the amount of the strain depends upon the external 
conditions of restraint or, in the case of a work-hardening material, on the 
vate of work-hardening. 

Usually, however, we have non-uniform stress distributions, and the amount 
of deformation depends upon the stress pattern as well as the boundary con- 
ditions. 

Der SAINT-VENANT (1870) [13] stated the law that the directions of the 
increments of the principal strains coincide with the directions of the principal 
stresses. In this statement DE SAINT-VENANT recognized that linear relation 
between stress and strain does not exist in plasticity as it does in elasticity. 
With plastic deformations, which may be very large, the general displacements 
of the material may be such that the directions of principal stress throughout 
the plastic mass will, and generally do, change. So we must follow the increments 
of plastic strain which result from a given stress condition. 

How big are the increments we are talking about? In theory, they are infi- | 
nitesimal since theoretically there is no difficulty in passing through an infinity 
of stress states. In practical cases, the increments will be anything from 0-1 per 
cent strain to 1 or even 3 per cent strain. A crude criterion is to say that the 
increments can be as large as can be considered not to change the stress distri- 
bution to such an extent that the principal stresses change direction. 

Lévy [14] first proposed that, over these strain increments, the strain is 
proportional to the stress. In terms of the deviatoric stresses, the relations are 


ved P. 
Stag) Sey ee (6) 
Se Bs Py 


These relations were restated by von MisEs many years later [15]. They are 
usually called the Lévy-Mises equations. 


Where both the elastic and plastic strains have to be allowed for—and strictly 
they should be in all cases—Reuss [15] put forward the following equations 


bez, = de8, + de?, = Pe +2 p'.... (three equations) , (7) 


Ol gy = ore +Abyy ++ (three equations) . (8) . 


a 
J 
‘! 


| Vol. V, 1954 Theory of Plasticity 7 


In all these equations 6 before the strain indicates that a strain increment 
must be taken, not the total strain. The factor / is an artifice which represents 
_the plastic strain increment as proportional to the deviatoric stress: it is not, 
_however, a constant as in elasticity but varies not only from moment to moment 
but from point to point. 

Using the Reuss or the Lévy-Mises equations depends upon whether the 
plastic strains are comparatively large or small, or whether the elastic material 
acts as a restraint to the plastic region. Where the elastic strains are ignored in 
the latter, there is a likely discontinuity since at the plastic-elastic interface 
a sudden change of strains will occur because the lateral contraction ratio 
will change over from Poisson’s ratio of about 0-3 in the elastic range to a 
value of 0-5 in the plastic range. 


Condition of Constant Volume for Change of Shape 


This last statement means that we assume constancy of volume in the 
plastic range 
det, + def, + der, =0. (9) 


Accurate measurements on materials heavily deformed plastically show this 
assumption to be very closely supported. Apart from material consolidation 
owing to internal cracks, blowholes, etc., materials do not change volume by 
more than a very small fraction of one per cent. 


Work-Hardening 


We can determine for pure shear or simple tension, the work-hardening 
curve of the material, but we must now ask the question ‘with a complex stress 
state, what is the new value of & after a given set of increments of the plastic 

strains?’ We assume that the hardening takes place isotropically and the usual 
answer is that the new & depends upon the work done, i.e. 


Jy = FW?) = F( | p* de"), (10) 
where ~* = equivalent stress, 
de* = equivalent strain increment. 


These work functions can be evaluated from the stress-strain curve in 


simple tension or in torsion [16]. im 
An alternative and usually more useful basis is the total amount of strain, 1.e. 


Je = H([ dex). (10a) 


Either basis appears to give satisfactory results. 
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Solution of the General Problem 


In general, we have nine unknowns ; six stress components and three displace- 
ment increments. We have nine relationships; three equilibrium equations, 
constancy of volume, yield criterion and four (independent) stress-strain in- 
crement relations (Lévy-Misgs or Reuss). Although there are six of these, two 
are inherent in the other relations. Alternatively we can use the statement that 
the principal stress and strain increment axes coincide and the three stress- 
strain increment relations to make four as before. 


Elastic-Plastic Problems 


The complete solution of a general problem in plasticity involves a calcula- 
tion of the stress and deformation in both the elastic and plastic regions. 

The solutions in the elastic and plastic regions are interrelated by certain 
continuity conditions in the stresses and displacements which must be satisfied 
along the plastic-elastic boundary. The boundary itself is usually one of the 
unknowns and is usually of such an awkward shape that even the stress distri- 
bution in the elastic region can only be obtained by laborious numerical methods. 
The complete and accurate solution of the plastic problem is practicable in 
relatively few cases. 

HOHENEMSER [17] attempted an experimental confirmation of the Reuss 
equations. The results were not conclusive, but the agreement was probably as 
good as could be expected in view of the uncertainties of this sort of experiment: 
hysteresis, the rounding of the yield point, anisotropy, etc. Such secondary 
factors tend to prevent these tests being critical when the strains are small. More 
recent tests by MorRISON and SHEPHERD [18] were in substantial agreement 
with the Reuss equations. 


Plastic Bending under Tension 


We shall first consider a problem in which the plastic and elastic regions 
both exist. The simple punch and die assembly, Figure 3, starts with a flat 
blank and as the punch penetrates the die, so the blank is drawn over the die 
radius into a cylindrical form around the punch, until all the radial blank has 
passed over the die radius into the die throat. 

SwirT [19] has found that the amount of thinning of the cupwall is much 
greater than would be expected from simple drawing theory, and has shown that 
this arises from bending the plate plastically (at A and B) under a tensile load 
in plane strain. (Plane, because the sudden change from a flat plate to the die 


radius and vice versa occurs with zero circumferential strain). Consider a beam 
undergoing plastic bending (Figure 4). 
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We cannot assume that the neutral axis NN coincides with the mid-fibre CC 


because there is an out-of-balance tension over the section corresponding to the 
radial pull P. 


P= | 4, dy. 


‘The distance yy will depend upon the extent of the plastic and elastic regions. 


Strain distribution 


—2k 
Stress Wstributron 


Figure 4 


Over the elastic region y = y, to > 
bp=a+by (straight line law). 


At the limits of the plastic zone ~, = 2 k, the yield stress in plane strain. 
We shall assume that the strains are small and that the yield stress can be con- 
sidered constant. 

Then ; 

ee, Vor Va Wee uaa ° 
ES oi ee a 


whence 
P=4k(d— yy). (11) 
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The limits of the elastic range are fixed by 


Zk Yea— In _ Ya— YN (12) 
E e R-+ yn 


where g = radius of curvature of neutral axis. 
Thus if the tensile load per unit width P and the die radius R are known, 
yy, y, und y, are known. If P is large and RF small, as is usual in deep drawing, 


the elastic region becomes so small that the stress distribution is essentially as 
in Figure 5. Then, in plane plastic strain 


6€, + dé, = 0 [equation (9)| 
and 
‘ : 2d— 4kda+ FP. 
de = thickness strain = : ees {ECL eR SS 
on the tension side, and 
YN 4 d—P 


Ugo ahd ae 


in the compression side, that is, a total thinning of 


_4#kd . 
4hd+4kR—P- 


(13) 


Plane Strain and the Slip-Line Field Theory 


By far the most frequent and important applications of plasticity theory 
involve relatively large plastic strains, and in this case the elastic strains are 
negligible both in magnitude and in their influence on the stress and strain 
patterns. 

So far, few general problems have been solved, even approximately, but for 
plane strain problems the theory of the slip-line field allows many cases to be 
solved. Plane strain conditions hold in many metal working processes: sheet, 
strip, and wide bar rolling ; bar drawing where the width is large compared to the 
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depth; and certain aspects of deep drawing—to quote a few examples. In all 
these cases, the plastic strains are large. A rigid plastic material is assumed and 
the Lévy-Mises equations are therefore used. Otherwise, in plane strain the 
most complete solution of plastic flow problems have been realized. 
The conditions are:— 
(1) Flow everywhere parallel to a given plane. 
(2) Motion independent of co-ordinate perpendicular to the given plane. 
Using co-ordinates x and y, these conditions mean that we have five un- 
‘known quantities. 
Stresses Dye, Py y, Pzy- Displacements u,, vy, and there are five equations :— 
Equilibrium Equations (2) 


be ey (14) 
OP, Op, 
ey ee (15) 
Criterion of yield (1) 
(Paw — Pyy)® + 4 Diy = 4 R?. (16) 
Constant volume (1) 
: Ou, ov 
es —— —lyy 1. Wve = ae 7 (17) 
Lévy-Mises equations (1) 
OUy, Vy 
et tye Cos at 
fe Pit Jin, Se re) 


It is possible to go into considerable detail with these equations, to derive 
various theorems and conditions which must apply. These have been dealt with 
exhaustively by Hirt [20] and by PRAGER and Honce [21]. There are some- 
times enough boundary conditions in the stresses to solve the problem without 
reference to the displacements. These are called ‘statically determined’ prob- 
lems, but there are not many of them. Usually the stress and displacement 
equations have to be worked together. Certain conditions about the stresses are 
deducible, however. Thus the stress normal to the plastic-rigid interface and the 
shearing stress on the interface must be continuous for equilibrium across the 
interface. Because the yield criterion is satisfied on both sides of the boundary 
(the rigid side is just on the point of yield) the direct stress in the tangential 
direction may or may not be continuous. This follows from the quadratic form of 
the yield criterion. 

Thus, if £,,, is the stress normal to the rigid-plastic boundary and the yield 
stress k is known, then: 


(Pez ~~ Pov) = 2) Be? Pav Le. PDyy = Pew a2 VR > Pay c 


| 
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In the limit, two values of p,, differing by + 2% can exist on either side 
of the boundary, when ,, becomes either the greatest or smallest principal 
stress p, or f, (Figure 6). The interpretation is to be seen in the example of | 
plastic bending under tension already discussed, where it was shown that, | 
when the radius of curvature is small compared with the thickness of the mate- | 


; 


Figure 6 


rial, the elastic region reduces to an imperceptibly narrow layer between two 
plastic regions, one side being tension, the other compression. PRAGER calls this 
‘the last vestige of an elastic region’, and since we are neglecting the elastic 
stresses, we can, in the extreme case, have a discontinuity of 4 k. 

The diagram Figure 6 gives a clue to the representation of the plastic con- 
dition at any point in plane strain. 

The yield criterion is 


un Ta buy)” ge Day = 4 


and this is the equation of the Mohr circle, radius k, centre at distance 
(Paw + pyy)/2 from the origin. The points J and JT then represent the yield con- 
dition and if we choose our element so that its sides are in the directions J and JJ , 
then the element corresponds to the directions of maximum shear stress, the 
normal stress being the hydrostatic stress (Figure 7). The sides are therefore 
lines of shear (slip lines) and all elements oriented in this way represent shear 
elements with no extension of their sides, which are of course orthogonal. 
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Velocity Components 


When dealing with slip-line fields, it is more convenient to consider velocity 
components, that is, the displacements wu, and v, are to be considered as the 
displacements in unit time, corresponding to the steady displacement of the 
boundary in the same time interval. Since none of our quantities are time- 
dependent (see assumptions page 2) the time rate of strain need not be considered 
as influencing the results. 


Figure 7 


Similar arguments can be applied to the velocity components for conti- 
nuity as to the stresses, but H1Lt has dealt with this matter in detail. It is only 
necessary to point out that the directions of maximum shear rate coincide with 
the directions of maximum shear stress. Such boundaries (and they form the 
slip lines of our plastic field) may represent discontinuities in the deformation 
pattern, just as they may in the tangential stress component. 


Curvilinear Co-Ordinates 


In plane strain, therefore, the slip lines will form a set of orthogonal curvi- 
linear co-ordinates throughout the plastic mass and the solution of any pro- 
blem is the determination of the equations for a curvilinear element. The equa- 
tions of equilibrium require restating in terms of these co-ordinates rather than 
the Cartesian co-ordinates of equations (14) and (15). 

For a very small element the equations of equilibrium become 


op a. 
59 2 k=0 (19) 
and on integration give 


p + 2k 0 =constant along a first shear line, (20) 


p — 2k 6 = constant along a second shear line (21) 
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where 0 is the angle made by the shear line to some reference direction 


(HENCKY [22)). ; ye 
The value of the constant is generally different from one slip line to the 


next, but along any slip line, if we move through an angle 6@ then 
p+2k0=const, p’+2k (6+ 60) = const 
whence 
Pp =p-—2k60, 
or, on the second slip line, 
p=p4+2k00. 


Velocity Relations 


Considerations of continuity and the conversion of the Cartesian compo- 


nents of velocity to curvilinear co-ordinates leads to the Geiringer equations. 
du+vdd=0, (22) 
dv —u dd =0 (23) 


where # and v are the velocities in the slip-line directions. 


General Theorems and Characteristics of Slip-Line Fields 


PRAGER and HopGE [21] have described many of the special properties of 
slip-line fields and it is unnecessary to go over the ground here. Usemay be made of 
any special characteristics which can be discerned when solving particular prob- 
lems, but reference must be made to a graphical method which PRAGER?) [43] 
has lately put forward for solving problems and which greatly simplifies the 
work. In this method the stress plane, physical plane (slip-line field), and velocity 
plane are constructed together. A few examples will illustrate this. 


Example 1—Extrusion of a flat sheet (Reduction ratio 3:1) 


The slip-line field is shown in Figure 8 (see PRAGER and HopceE, page 182) 
for a half section through the die (the other half is of course a mirror image). 
The die is taken as perfectly smooth, so that there is no friction at the die wall 
and the slip lines meet it at 45°. Material to the left of the limiting slip line 7 to 8 
moves at uniform speed “, that beyond the line 7 to 0 moves at uniform speed 


1) Originally presented in Sweden, 1952 [43], but first given in detail by PRAGeErR at the 
Summer School on Plasticity at Imperial College, England, July 1953. The author wishes to pay 


a personal tribute to the considerable influence which PraGeR has exercised in every phase of 
Plasticity. 
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Physical plane 
(slip lines) 


Stress Plane Fs 


Absolute 3 a f 
velocity “ 


Figure 8 


3 u, to the right because the reduction ratio is 3:1. At both slip lines there is a 
discontinuity, and considering any element A arriving at velocity u, at the 


entry to the plastic zone, there is a velocity discontinuity as it crosses the line. 
The Geiringer equation (22) gives 
dv —udd=0. 


Relative to the point A on the left of the slip line, the element has velocity 
components # in the normal direction and v in the tangential. From continuity 


considerations, w’ on the right of the line must equal uw (Figure 9), i.e. 


dv=u0dd, dv’ =u' d0=u 00. 
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Therefore, as we muve along the slip line through any angle 06 


dv = ov’ 
1.e. 
v —v’ =a constant along a slip line, 
= the discontinuity at the slip line. 


If we can find the velocity discontinuity at any point, then we know it all 
along the slip line. Referring to Figure 8, at point 8 an element arrives at velo- 


Velocity discontinuity 
=e, (slip line) 


v'-v-Amount of 
discontinuity 


Path of element 


Figure 9 


city uw, and must turn sharply through a right angle: it must move downwards 
at the same velocity (constancy of volume) because the velocity in the x-direc- 
tion is zero the moment it passes across the limiting slip line, because of the die 
wall. We start to construct the velocity diagram. From pole point O, length 
OA = uy the velocity of approach, length 0 to 7’, 8’ is the vertical velocity of 
any point on the right of the slip line 7 to 8 and the length A to 7’, 8’ equal to 
V2 u, is the discontinuity along the slip line. Within the triangle C, the material 
moves down uniformly as a rigid zone until the line BC is reached. The 
region B is a centred fan, with a discontinuity of value /2 u, all along the arc 
AB. Any element, 5 say, has velocity u, horizontally, /2 uw, tangentially to the 
limiting slip line, and the arc points 7 to 7 are represented on the arc 7’ to 7’ 
_on the velocity plane. At 7, or at any point along the limiting line BD, there 

is a further discontinuity, and to complete the velocity diagram so that we reach 
point D, the discontinuity at the line is V2 u,. 


Streamlines 


The absolute velocity of any point is represented by the vector joining the 
point to the pole; e.g. a point immediately to the right of the limiting slip line 
at 5 has velocity 0 to 5’. The streamline from A can then be drawn in, by draw- 
ing short chords parallel to the absolute velocity at each successive point. The 
streamline #, g, s is shown. 
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Example 2—Compression between smooth platens 


If a slab of an ideal rigid-plastic material is compressed between smooth, 
tigid platens (Figure 10a) and the width of the tool, 5, is equal to the slab 
thickness ¢, then the slip-line field is a simple criss-cross of shear, all deformation 
being concentrated on the two 45° lines of discontinuity, and if k is the shear yield 
Stress in plane strain, then it can easily be shown that the mean pressure on the 


ASB CORrae 


f 
Die breadth, b 


Figure 10a. Solution for t ~ 0-7 b. 


die ~, is 2k. As the thickness ¢ becomes progressively greater, so the slip-line 
field changes until the slab is so thick that the field from the top die no longer 
influences that of the bottom one, and it can be shown [23] that the mean die 
pressure in the limiting case is 


po= 2k (1+ 5) =2571 x 2k. (24) 


This limit is reached according to Hii when ¢/b > 8-74. 

_ If we start again at ¢/b = 1 and imagine the dies to approach each other 
steadily, the slip-line field will again go through a gradually changing pattern 
until we reach a ratio ¢/b = 1/2. Here the slip-line pattern becomes a double 
criss-cross of 45° lines, and this pattern will be repeated at ¢/b = 1/3, 1/4, etc. 

GREEN [24] has examined the slip-line fields for a number of points between 
(/b = 1 and ¢/b = 1/2, and from the stress plane the mean pressure py on the die 
can be calculated. Figures 10b-d shows the physical (slip-line) plane, the stress 
plane, and the velocity plane for the case where ¢/b ~ 1/2. 

Having calculated a number of such cases, GREEN showed that, as the dies 
srogressively approach, from condition ¢ = ) to t = b/2, the mean pressure rises 
0 a maximum at ¢ = (1///2) b some 4 per cent higher than the yield stress 2 k 
n plane compression, falling again to 2 # as it reaches ¢ = 2 b (Figure 11). 
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Figure 10b. Solution for t ~ 0-7 b. 
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Figure 10d. 
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Figure 11 


Experimental Investigation of Slip-Line Field Solutions 


| It may be asked, how far do we know from experimental measurements that 
these slip-line field results are correct? 
There are hardly any experiments which are conclusive because most in- 
-vestigators have used small deformations and have tried to observe the first 
yielding of material, which is so much dependent upon the elastic behaviour 
and the transition from the elastic to the plastic region. 
_ However, for several years the author has been working on plane plastic 
‘compression, using highly polished and lubricated dies having a width which 

does not overlap the material (Figure 10a). The unconfined material on each 
‘side of the dies acts as a very rigid lateral restraint on the material between the 
dies and even for quite narrow strip, the condition corresponds remarkably 
closely to plane strain [25]. 

This test can be made to give conditions remarkably close to those of the 
theory and indeed, while GREEN was working on the theoretical problem des- 
cribed above, experiments were being carried out at Imperial College quite 
independently over the same range of geometrical conditions. If a techno- 
logical material is to behave as an ideal one, then certain characteristics must 


be aimed at: — 
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(1) There must be no anisotropy. 

If we use cold rolled high conductivity copper, the strain has taken place 
in the same sense that the further strain is to occur in plane compression and 
the anisotropy is of little importance. Tests have shown that the materia] has 
rotational isotropy in the vertical plane. 

(2) The material must have a sharp yield point. 

If we suppose the material cold-worked to point F (Figure 12) and then 
tested, say in plane compression, a curve FDBE would be obtained. There are: 


Nield-sfress curve 


0 e Strain 
Figure 12 


several points which might be used to give a value of the first plastic yield: 
D is the point where plastic deformation of at least part of the material has first 
started. B is the point where the curve from F becomes concurrent with the 
‘yield-stress curve’, that is, the curve that is the envelope of all stress-strain 
curves such as F’. It is difficult to determine either of these points and neither 
can be truly called the yield stress at the point F. A proof stress c is hardly 
appropriate in plasticity problems. 

The best method of determining fy is to extrapolate the curve EB to A. A is 
the point on the ‘yield stress curve’ corresponding to zero deformation. It is 
a clearly defined point and has the advantage of corresponding to the theoretical 
condition of no deformation, eliminating the elastic strains. 

(3) There should be no work-hardening. 

If high conductivity copper is work-hardened by cold rolling, the yield-stress 
curve beyond about 30 per cent reduction in thickness, becomes nearly straight 
and shows little work-hardening. If we can extrapolate to zero deformation 
each time, the work-hardening is less important, but if we attempt to measure pp 
the mean pressure, using one pair of dies of width b, then ¢ must vary and each 
stress state will be super-imposed on the previous one, and the effects of flow 
cannot be separated from those of yield, because in effect we are moving the 
point F up the yield-stress curve. In these experiments } was changed and ¢ 
held constant, so that the yield stress of the material remained constant for 


all the tests, and the dies were simply moved to a new point on the same strip 
for each die width. 
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(4) There must be no friction. 

This condition cannot be exactly met, but this requirement and requirement 
(2) were very closely approximated by an incremental loading method, which 
allowed the extrapolation to be made. 

The extrapolation was carried out in the following way: — 

The strip and dies were measured by micrometer, a pencil line across the 
strip being used to indicate the position at which the indentation was to be 
made. The strip was then lubricated at this point by the application of a soft 
graphite pencil, with a touch of mineral oil, over the whole area. This method 
of lubrication has been found to give more successful results than pure graphite 
flakes and has the advantage of being very convenient and simple. After lubri- 
cation, the strip was inserted between the dies with the edges at right angles 
to the width of the dies and centrally placed with respect to the loading axis. 
Load was then applied to the subpress until it reached a value necessary to pro- 
duce a small deformation. At this value the load was held steady for 15 to 30s, 
then released and the strip removed from the subpress, the mean thickness of 
the strip at the indentation being then measured. The strip was lubricated again, 
replaced carefully between the dies to register in the indentation already made, 
and a slightly greater load was applied. This ‘incremental’ loading process was 
repeated several times, in the same indentation, until enough of the curve had 
been plotted to allow the extrapolation to be made. The load increments chosen 
in the tests were large enough only to ensure that the load-compression or ‘yield 
stress’ curve BE was reached each time. 

This ‘incremental’ loading method overcomes the difficulty of measuring the 
strip thickness, because the strip is removed between load increments. Another 
advantage of the method using one indentation is that the strip can be lubricated 
between each increment and, since these increments of compression are small, 
the friction cannot develop its full effect and is kept at a very low value indeed. 
The conditions, therefore, come very close to those specified in theory. 

If we plot p)/2 & against the ¢/b ratio, Figure 11 is obtained, where GREEN’S 
theoretical result is compared with the experimental. The agreement is remark- 
ably close. The experiments extended to ¢/b = 1/3, and GREEN has calculated 
the upper and lower bounds [26] of the oe fields and our results over the 
next cycle also fit nicely. 

In the range where ¢/b > 1—that is, in the range where Hixv and others 
have similarly calculated the fields—the agreement is equally good (Figure 13). 


The Rolling of Metals?) 


_ The problem that all theories of rolling aim to solve is this:—given that it 
is desired to reduce the thickness of a piece of material by a certain amount, 


1) In the following sections, compressive stresses are reckoned positive. 
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Figure 13 
Curves showing the variation of p,/2 k with t/b (t/b >1). 


what is the nature of the deformation and what roll force and torque need to be 
applied to the rolls? The information needed to solve this problem is: the 
initial and final thickness of the strip, the yield-stress curve of the material, 
the radius and elasticity of the roll, the coefficient of friction, and any exter- 
nally applied forces such as strip tensions. 

In the theory of plasticity, the assumption is made that the material’s 
resistance to deformation is independent of the rate of strain. This is the con- 
dition during cold working, that is, where the material does not recrystallize 
faster than it is deformed. In hot working, however, the yield stress of many 
materials is dependent on rate of deformation and another factor enters. 

Consideration will, therefore, be confined to cold rolling, in which the yield 
stress only depends upon the amount of strain, and to strip or sheet rolling, 
where the width is enough to ensure plane strain conditions. The problem in 
some respects is very much the same as that of plane compression between 
parallel dies, but the ratio of thickness to die breath is much smaller and the 
dies are slightly inclined to each other. The theory of the slip-line field for 


inclined dies has been evaluated but it has not been applied to the case of 
rolling because: 
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(1) The curvature of the dies will have an unknown influence. 

(2) Although in cold rolling the coefficient of friction between roll and strip is 
small, it has a great influence on the process: indeed the process is not pos- 
sible without it. 

An attempt was made by OrowaNn [27] to reach a more accurate solution 


- of hot and cold rolling, but for cold rolling, the assumptions made by. von 


AV PEE Chee 


KArmAN [28] are adequate, although the latter’s solution is not satisfactory for 
hot rolling. 

When strip is rolled between smooth, lubricated rolls, the strip is in equi- 
librium under the effect of any external tensions, the roll torque being such as 
to cause the vertical reduction in thickness. The reduction in height causes a 


x 


_—_—_———m— 


Figure 14 


corresponding elongation, so that the speed of the strip leaving the rolls is 
greater than that entering them. A balance is set up, the rolls rotating at a speed 
which is less than that of the strip leaving and greater than that of the strip 
entering. So there is relative movement between the strip and roll surface, and 
since there are high radial pressures causing the deformation, there are frictional 
forces developed, and these act in opposite directions, depending upon whether 
we are on the entry or exit sides. At some intermediate point, the ‘neutral point’, 
the roll and strip must have the same surface speed and the direction of the 
tangential force changes over. This means that the radial force must be increas- 
ing as we move inwards from entry to the neutral point and likewise from exit 
giving a roof type pressure distribution. This is called the ‘friction hill’. The 
horizontal component of the frictional forces and of the outward component 
of the normal force must just balance the resultant external tension force. 

Von KArMAN made one main assumption: 

Plane vertical sections remain plane, i.e. there is slipping friction throughout 
the contact arc and the vertical and horizontal stresses are principal stresses. 

There were also two approximations: 
(a) The vertical component of the roll pressure, g, can be taken as approxima- 

tely equal to the roll pressure S. 


(6) The arc of contact is circular. 
These conditions are illustrated in Figure 14. The vertical AB remains 


straight and vertical as it passes through the roll gap, and any friction there is 
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at the roll surface has no effect in modifying the stress and strain patterns. The — 


deformation is homogeneous; naturally in cold rolling the angle of contact is 
not nearly as large as shown and under these circumstances the approximation 


Figure 15 


is reasonable. VoN KARMAN’s theory, therefore, really only applies to cold 
rolling. In hot rolling, there is high friction and the material is sheared over at 
the surface, as shown exaggerated in Figure 15. 
If the angle of contact is small—and in cold rolling it is usually less than 5°— 
it is obviously reasonable to take the normal and vertical components as equal. 
The vertical element (Figure 16) will be under the external stresses, S, uw S 
at the surfaces and f, p + 6p on the vertical (internal) faces. We assume that 


Elastic entry Pla 3 ’ _ | Liasfic ey, 


rolling =p 


Figure 16 


the arc of contact remains circular, radius R’ ; thus for an element at angle Q, 
width R’ dy, the horizontal equilibrium requires 


(p + df) (4 + 6h) —ph=2SR dpsing+2uS R dp cosp 


| 
: 


a en ae, oe 
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1.e. 
PRs S (sing - uw cosq) . (25) 


This is VON KARMAN’s equation. 
By the Mises criterion © 


s—p=2k (assumption 2) 


where 2 kis the yield stress in plane strain = (yield stress in simple tension) x 2/V3. 


Solutions of von Kdrman’s Equation 


All the argument about rolling theory really reduces to solving equation (25) 
and the differences are in the approximations or further assumptions made. 
Many solutions so mutilate this equation that they are no more than empirical 
results without any real scientific basis. 

OROWAN [29], without approximation (a) showed that a stepwise numerical 
solution was possible, allowing for k varying through the roll gap and was the 
first to give the correct solution of the equation (called the homogeneous gra- 
phical method). It is, however, rather laborious in use and it can be shown that 
certain further approximations can be made which simplify calculations, yet do 
not sensibly reduce the accuracy of the results [30]. 

- We make the usual assumption that the coefficient of friction, wu, is constant 
over the arc of contact. 

For angles less than 5° 

sng, cosp=i, 


except where (1 — cos @) occurs multiplied by a large factor when 


2 
1 —cosp = -. 
For example 
h=hg+2R’' (1— cosy) =h,+ R'g?. 


The roll gap is divided into three zones (Figure 16): 
(1) Elastic compression at entry. 
(2) Plastic compression. 
(3) Elastic recovery at exit. 
We shall first ignore (1) and (3) 
The geometry of the roll gap then becomes: 
length of contact are 
L=R! (hy — hy) , (26) 
angle of contact 
fia ees Be (27) 
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Then from the differential equation (25), i.e. 


AUD (2 1) +hh ge (S)=2R's tH). (28) 


The boundary conditions are given by the front and back tensions. At entry 
P= h=hk, \—(4)=h, 5=%—4 (total tension = hy ty) . 


At exit 
p=0, h=h, s,=hky—t, (total tension = hp t) . 


The solution builds up from each end and the curves meet at the neutral point py. 
Consider first the case of no tensions. 
We assume that the product 4k does not change very rapidly: for most 
metals this is so; 1.e. 
k(hk aa [oS! 
“ig SPP ale) (29) 
Therefore, 
a fS 
Ten _ 2R’ (pth) 
Sk Ie + R’ g? 
On integration, therefore, 


sh gad 
S= Cie (30) 
where 
Re we 
. H=2 |/ ie tan (|/ i °) ; (31) 
At exit 
So = %s, at entry S;= x, - (32) 
On exit side 
kh 
pe 6 (33) 
ron hae a 
aT aed (34) 


where H, = value of H at p = gy. 

The above approximation (29) does not hold good with annealed strip and 
high back tension. But the equation can easily be modified for the case where 
tensions are applied to the strip [31]. We imagine two passes, one in which no 
ee is applied and another identical with the first but with applied tension ¢, 
and fp. 

Call the roll pressure S, in the no-tension case, then from the differential 
equation 


Sq dk ; 
dp = 2uR So—h | —2R’ kp =0 


ee 


SA ae Oe Sa oe ee 
. \ ™ 
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and for the other case 

dS . dk ; 
subtract and integrate 


d 
h Gy (So— S) = 2" R’ (Sp— S) =0 


or 
SiS Sal Seg (35) 
At entry 
(So — S,={f ed (ace ‘ij }= fle eae Sie 
At exit 


(So — S)p={ ky — (h, —#,) }== 4, e707 1. 


Inserting in (35) we get B* and therefore 


St= G- = t) et”, (36) 
s-= G- Spent (37) 


These equations are satisfactory over the whole range of tensions and pass 
reductions. 


Total Roll Force 


For the total roll force, we integrate the area under the curve of S. 
The roll force per unit width P is 


Pi Pi ¢ 
p=/SRdp—R' | sap. (38) 
0 0 


This must be done numerically by Smmpson’s rule or other area formula, 
because S depends on k, which varies with degree of reduction. 


Roll Flattening 


The force P acting on a small area of the roll surface causes a flattening of 
the roll. . 

Hircucocx [32] proposed that if the roll pressure distribution could be con- 
sidered as approximately elliptical—and this is not far wrong—then it can be 
shown from the Hertzian theory of contact that the radius of curvature of the 
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deformed roll is constant over the arc and equal to 
’ cp 
R'=R(1+“F) . (39) 


where R = roll radius, 
P = roll force per unit width, 
¢ =a constant = 3-34 x 10-4 for steel rolls, but depends upon the 
elastic properties of roll material, 
6 = the pass reduction = (A, — fy). 


Roll Torque 
OrowaNn [29] showed that the torque per roll, per unit width 


Pr Pn 
c= | s-ay—['s-a). (40) 
0 


Pn 

This relation is the difference of two nearly equal amounts, and with an 
approximate method of solution, is not accurate. It was shown [30] that an 
expression insensitive to slight inaccuracies in calculating the pressure distri- 
bution could be derived from the equilibrium of the horizontal forces: 

Pi Gn Fr 
(thy — ty) —2 fw S- Rdg +2 | St R'dp +2] S@R' dp =0 
n 0 0 


whence, using (4), therefore 
Pr 


G=RR'| Spdp+ > h-hh). (41) 


0 


The Elastic Arcs of Contact 


Although it was realized at an early date in the development of rolling 
theories that the roll deformed elastically, no allowance was made for the effect 
of the elastic strains in the strip. For large passes on fairly thick strip the 
elastic arcs are comparatively small and the error is unimportant. 

For medium and small passes the elastic recovery at exit has a sizeable 
influence and in very small passes—4 per cent reduction or even less—even the 
elastic compression at entry may be significant. 

The total roll force 


Diss Poe dai Eee (42) 
Considering only the elastic recovery zone (Figure 16), (e. P =P; 4B) 
at the end of the plastic zone and beginning of the elastic zone the tension is tes 


and the minimum height /,,. At true exit f, and hy respectively. We take as a 
first approximation that the elastic recovery zone is rectangular, of height h, 


oe RL eS Oa Oey ee oa elec i saa. ie 
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_ and length a. Since the coefficient of friction is usually not much greater than 


0-06 in commercial cold rolling, the frictional force over the surface of the 


_ elastic arc is small and we neglect it. » 


We do not, however, neglect the effect of the frictional force in calculating 
the effective strip tension ¢,, acting on the plastic zone neither do we ignore it 


_ when calculating the roll torque. 


Using these assumptions and approximations it is possible to calculate the 


_ overall elastic contribution to the roll force in terms of the known displacement 


of the strip surface in the elastic region [33]. The length of the elastic recovery 
arc, a, is given by: 


Ine ae 
ea ET (1 — »?) (he, — ts) 


and the contribution to the roll force P,, 
P., = (a constant) x V ho Rilke iy (43) 


where E = elastic modulus for the strip and v is Porsson’s ratio. 
In terms of the externally applied tension ¢, 


2"P, 
ly = loo + i . (44) 


The roll torque is modified by the amount 
Geo = —- R Pee F 
The radius of the flattened arc of contact, R’, is slightly modified [33] 


Ge, 
R= Ri 1 ee —— \, 
{14 (Vea 8 + 8, + os)? 
fs) = hy = Ng, 
ks (1 — v2 
0, = h, — hy = yo ie me (ha — tes) 


(h, = thickness with tension still applied), 


6,= ena, (Ig ty — hy 4.) = recovery at exit when tension removed. 


Approximate Methods of Calculating Roll Force and Torque 


There are so many calculating methods available that the engineer may 
very well feel uncertain as to whether one method is better than another but 
there are several questions which may be asked in deciding on a reliable 
method to use. 

For general use, the real basis is to use only the yield-stress curve of the 
material to be rolled and the coefficient of friction between the roll and material. 
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The solution of von KARMAN’S equation already described allows for the varia- 
tion of the yield stress through the roll gap, and necessarily involves fairly 
lengthy calculations. A simplification for practical calculations would be pos- 
sible with a constant yield stress because equations (36) and (37) are such that 
the neutral point can be predicted and the shape of the pressure distribution 
depends only upon the shape of the yield-stress curve. 

The early theories of rolling all assumed a constant yield stress. That the 
predicted results were so much at variance with each other and with the miscel- 
laneous experimental results then available is not surprising, especially as these 
early theories were developed for hot rolling, where the assumptions made do 
not hold. 

The first requirement for a practical calculation method is to make equa- 
tions (36) to (41) independent of k. For the case of rolling without tension, CooK 
and LaRrkKE [34] showed that the simple Ekelund formula 


rar 16u/R’ d—1-26 
P= AVR’ (1+ ss (45) 


gave results in satisfactory agreement with OROWAN’s theory, except for a 
constant varying from about 0-98 to 1:09 depending upon the material. 

If we plot per cent reduction against position in the roll gap, it will be found 
that most of the deformation has occured in the first 40 per cent of the contact 
length, so the average yield stress is not the mean effective. Cook and LARKE’S 
correcting factor was found [35] to be unnecessary if the true mean yield stress 
was taken 

f f 
[R dp [dg 
k=* =< (46) 


= Pi Pi 


This ‘bunching’ of the deformation at the entry side makes it possible to use 
a mean yield stress, without much loss of accuracy. 

With a constant mean yield stress equations (38) and (41) can be numerically 
integrated once and for all over the practical range of cold rolling conditions [36] 


and graphically presented as a roll force function in terms of three non-dimen- 
sional quantities: 


as ; 
Gs le (geometrical factor), 
2 
b= pos (tension factor), 


= —>-—— (reduction aactor., 
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The roll force per unit width is then: 


P, = (k—4,) VR’ 6 x p(a, 7) (47) 


p(a, 6, r) being the roll force function. 
A similar function for torque can be obtained from equation (41) 


G=R6(k—2) x g(a,r,d). 


Although a double interpolation is necessary, this can be done quite quickly 
in practice and a further refinement is to allow for the roll force and torque 
resulting from the elastic recovery [36]. The method is simple, but necessitates 
a trial and error calculation since the elastic arc causes a compressive force on 
the end of the plastic zone and this must be subtracted from the external 
tension. 


Comparison with Experimental Results 


Many thousands of separate rolling tests were made on the Sheffield Experi- 
mental Mill, using various pass reductions from 2 per cent up to 95 per cent 
in a single pass, with various types of materials, various lubricants, strip thick- 
ness, roll speed, and strip tensions. One of the most difficult factors to control 
is the friction in the pass, because most lubricants are very sensitive to small 
temperature changes and the roll force and torque change continuously during 
a pass on a long coil of strip [37]. It is of course particularly difficult to control 
when tensions are applied because tension variations cause the work done in 
the roll gap to change. The question arises—how accurately do our calculation 
methods and theories predict these results ? 

Two things influence the answer to this question. First, the effectiveness of 
the underlying theory in predicting the results (obviously our calculation 
method can be no better than the theory on which it is based) and second, the 
means we have for independently measuring yield stress and friction coefficient. 
If we accept that the yield-stress curve can be determined, then our only 
adjustable parameter is the coefficient of friction. We now have a little evidence 
about the friction likely to exist in cold rolling and it is possible, therefore, to 
test the calculation method against the experimental results [36]. 

About 150 separate passes made in the Sheffield Mill over a wide range of 
pass reductions and tensions were calculated using the method outlined above. 
Assuming the measured results to be correct, the percentage errors in the cal- 
culated roll force and torque were determined and it was found that 95 per cent 
of the results are within + 10 per cent of the experimental with passes varying 
from 4 to 74 per cent in a single pass. It is probable that the measured values 
may involve errors up to -+ 3 per cent, so that the calculation method appears 
to be very satisfactory. 
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Roll Flattening 


In all the foregoing calculations, Hircucocx’s formula was used to calculate 
the roll radius. This assumes that the pressure distribution can be taken as el- 
liptical, and while this is not a bad first approximation, it may be thought that 
the actual roll shape and the pressure distribution—being mutually influencing 
—will depart sensibly from so simple a picture. 

OrowaN pointed out that the pressure distribution clearly cannot be a 
peaked curve as calculated, The peak would cause a depression in the roll, and 
the made some tests which, while not inclusive, gave support to this view [27]. 

This effect is only found in very thin hard strip where the flattening is very 
great. BLAND [38] calculated the displacements along the arc of contact so that 
the plastic stresses and strains in the strip were consistent with those for the 
elastic roll and found that the theory requires also that the roll should show 
such an indentation. H1rcHcocK’s formula was shown to be very adequate, 
however, for most practical cases of cold rolling. 


Yield Stress Curves 


We have assumed in the foregoing that the yield stress curve is known, not 
only to some fundamental measure of strain but also to a basis of percentage 
reduction in the pass. It will be clear that the complex stress-strain distribution 
in the roll gap will cause a certain amount of work-hardening in excess of that 
represented by homogeneous compression from a height h, to h,, because there 
will be redundant shearing not indicated by the pass reduction. 

What we have to do is to determine from a simple mechanical test the 
stress-strain curve of the material, for example, in simple homogeneous tension 
or compression and interpret this in relation to the actual deformations in 
rolling. 

The usual stress-strain curve is found for simple homogeneous compression 
or tension: this can be done successfully up to extensions of about 40 per cent 
with ductile materials but this corresponds only to about 30 per cent reduction 
in thickness and is insufficient as a basis for rolling mill calculations. 

The method first investigated by Sacus [39] and later developed by Coox 
and LarkE [40] gives a means of establishing the yield stress in compression, 
for thick material, but is not applicable to thin strip material. Two methods 
were devised for determining the yield-stress curves for the investigations 
referred to above [35]. 

One uses the plane compression test described on page 17. The roll surfaces 
are very nearly parallel in cold rolling, and since the theories assume homo- 
geneous compression, the yield-stress curve in plane compression should be very 
close to that in cold rolling. In the test, a strip of material is compressed incre- 
mentally, lubricating each time, the successive increments giving progressive 
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reductions. The envelope of all the individual compressions forms the yield- 


stress curve, and since the friction is very small, a very near approach to the 
true plane compression curve is obtained. The ratio of strip thickness to die 
breadth should be kept within a ratio of 1:2 to 1:4 for the best results. 

The second method uses the tensile test, test pieces being cut in the length 


- direction of the strip from material progressively rolled down, say in 5 per cent 


passes. The true stress-true strain curve from each test is set up at the appro- 


priate ordinate of reduction, and the envelope of all the curves drawn in as the 


yield stress curve of the material. The basis of this method was fully described 
elsewhere [35]. 


Coefficient of Friction in Cold Rolling 


What do we know about the coefficient of friction between lubricated sur- 


_ faces pressed together with high pressure, and moved relatively to each other 


for small distances? The evidence is contradictory at present. The values 


_ derived from rolling tests as described above depend upon the theory used, the 


constancy of the rolling conditions, and whether the elastic recovery of the strip 


is allowed for. From many hundreds of passes which he has calculated at one 


time or another or by one method or another, the author has come to the 
conclusion that, if a coefficient of friction as normally conceived exists under 
rolling conditions, then the value is rather lower than has usually been assumed: 
probably between 0-05 and 0-07 for good lubrications and adequate cooling. 
For good lubricants, OROWAN and Los [41] found: 
(a) For mild steel on hardened steel, using palm oil w started at 0-025-0-035, 
rose gradually to 0-06 as the relative amount of movement increased. 


(6) For aluminium on hardened steel and slightly contaminated surfaces, yu 


started at about 0-08 and rose steadily with relative movements. 
For the same combination, but using Shell lubricant, w rose from 0-06 


B fo 0-1. 


The point OROWAN makes is that the linear movement is a controlling factor. 
These values are roughly confirmed by rolling tests, in that steel-on-steel 


seems to lead to lower coefficients than non-ferrous metals-on-steel. Moreover, 


by taking certain limiting rolling conditions, and allowing for elastic recovery, 


mw is found to range from 0-05 to 0-08. The relative movement in rolling is 
usually small corresponding to the lowest values given by OROWAN (movements 


less than 1 mm), but with large diameter rolls the relative movements might 


be more than this. Yet so far as we know, the coefficient is the same for all 
- geometrical similar systems, whereas such a phenomenon as described by 


OROWAN would involve a variable of the system. 

BowbDEN and Tasor’s [42] tests, using a special type of apparatus, gave 
coefficients for well lubricated surfaces from 0-07 to 0-10, but their test proce- 
dure, like OROWAN’S, involved sliders rubbing on surfaces and the conditions 
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are strictly not the same as in rolling. There is much surface smearing in rolling, 
but little transfer of metal to the roll surface. 

One method of finding the coefficient of friction in cold rolling tests, which 
is independent of any particular theory of rolling, is being used at Imperial 
College at present. It gives good, consistent results. 

An increasing back tension is applied to the strip until slipping occurs, when — 
the neutral point has moved to exit, i.e. y, = 0. Under this condition’ 


V1 Pi 
P=R [say and G=[RR'usdp 
6 0 


whence = G/(R P). Experimentally, it is necessary to measure and record con- 
tinously the roll force and torque in a suitable mill, in which a gradually increas- 
ing back tension can be applied. It is also necessary, strictly, to measure the 
forward slip (i.e. the amount by which the strip speed at exit exceeds the roll 
speed, as a fraction of roll speed) but in practice the value of G/{R P) peaks 
distinctly at zero forward slip so that no difficulty is experienced in observing it. 

This method has now been tried out for a great number of lubricants and 
rolling conditions and has proved very successful. It confirms that the coeffi- 
cient of friction lies between 0-04 and 0-08 for normal roll lubricants. The results 
are not yet published. 
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Hautkondensation an feingewellten Oberflachen 
bei Beriicksichtigung der Oberflachenspannungen’) 


Von ROMANO GREGORIG, Ziirich?) 


1. Vorbemerkung 


Es bestehen bekanntlich zwei verschiedene Arten von Kondensationen: 
Tropfenkondensation und Hautkondensation. Die Warmeiibergangszahlen bei 
Tropfenkondensation ergeben fiir tibliche Betriebszustande Werte, die ein gros- 
ses Vielfaches derjenigen bei Hautkondensation sein kénnen. Es ist somit klar, 
dass die Tropfenkondensation bevorzugt wird, die jedoch nur bei ganz beson- 
deren Zustainden der Oberflache, auf welcher die Kondensation stattfindet, auf- 
tritt. Diese kann mit Fliissigkeiten von geeigneter Oberflachenspannung derart 
behandelt werden, dass sich, wenigstens fiir eine kurze Zeit, Tropfenkonden- 
sation einstellt. Es ist jedoch bis heute noch kein Verfahren gefunden worden, 
durch welches man auf die Dauer zum Beispiel eines Jahres Tropfenkonden- 
sation garantieren koénnte. Deshalb ist man vorlaufig gezwungen, samtliche 
Warmeaustauschapparate, bei welchen Kondensation vorkommt, unter der 
Annahme von Hautkondensation zu dimensionieren. 

Es ist bekannt, dass bei einem Warmedurchgang nur der grésste Wider- 
stand massgebend ist. Daraus folgt, dass eine Verbesserung des Warmeiiber- 
ganges auf der Kondensationsseite nur dann besonders zweckmiassig ist, wenn 
der Widerstand der Warmestrémung auf dieser Seite grésser ist als auf der 
kithlenden Flache. Dies ist besonders bei vertikalen Vorwaérmern von Dampf- 
turbinenanlagen der Fall, wo die Warmeiibergangszahl auf der Kondensations- 
seite nur einen Bruchteil derjenigen auf der Warmwasserseite betragt. Bei 
solchen Vorwarmern lohnt sich die Verbesserung des Warmeiiberganges auf 
der Kondensationsseite. 


2. Wirkung einer feingewellten Oberflache auf den Vorgang 
der Hautkondensation 


Um den Mechanismus der laminaren Hautkondensation besser tiberblicken 


zu kénnen, sei hier die von NussELT im Jahre 1916 aufgestellte und noch heute 
giiltige Formel ins Auge gefasst. 


1) Erweiterte Fassung des am 3. Mai 1952 im Rahmen der Frihjahrstagung der Schweizeri- 
schen physikalischen Gesellschaft in Baden gehaltenen Kurzreferates. 
2) Forschungsabteilung der Escher Wyss AG. 
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Definiert man: 
a = Warmeiibergangszahl der Kondensation, 
y = Verdampfungswiarme, 
y/g = Dichte des Kondensats, 
A = Warmeleitzahl des Kondensats, 
m = dynamische Zahigkeit des Kondensats, 
@ = mittlere Temperaturdifferenz zwischen dem kondensierenden Dampf und 
der diesem zugekehrten Wandseite, 
H = Hohe der vertikalen ebenen Wand und 
g = Erdbeschleunigung, 


so ergibt sich nach der Nusseltschen Theorie: 


* fot ANGE) orl 
a = 3400 { “wie (1) 


Diese Gleichung gibt uns an, welche Gréssen geandert werden miissen und 
in welchem Sinne dies geschehen soll, um eine Vergrésserung von « zu er- 
reichen. Der erste Bruch in der geschweiften Klammer kann unter Beibehaltung 
des gleichen Dampfes nicht gedndert werden, weil er nur Stoffwerte enthalt. 
Die Temperaturdifferenz # kann insofern nicht variiert werden, als sie von den 
Gegebenheiten des Warmeprozesses vorgeschrieben ist. Es bleiben somit nur 
noch die zwei Gréssen g und H, iiber welche man frei verfiigen kann. Wenn 
jedoch der Abtransport des Kondensates nur unter der Wirkung der Erd- 
anziehung stattfinden soll, so fallt auch g als Variable ausser Betracht. 

Die Warmeiibergangszahl « kann nach Gleichung (1) gross gemacht werden, 
indem H klein gewahlt wird. Wird ein vertikales Rohr mit einem Draht be- 
stimmten Durchmessers spiralférmig umwunden, so reduziert sich die in die 
Gleichung einzusetzende Héhe H auf die Ganghdéhe der Spirale. Man erhalt 
eine Erhéhung der Warmeiibergangszahl annahernd proportional der vierten 
Wurzel des Verhiltnisses zwischen der gesamten Rohrhéhe und der Ganghéhe 
der Drahtspirale. Konstruktive Griinde haben die Anwendung von solchen 
Rohren als unzweckmiassig erscheinen lassen. 

Ein Ersatz der Erdanziehung mit durch Rotation hervorgerufenen Trag- 
heitskraften fiihrt zu praktisch unméglichen Ausfihrungen. 

Werden in einem Kondensator (in den Rohren das Kiihlmittel, aussen kon- 
densierender Dampf) nach Figur 1 vertikale, aussen gerillte Rohre vorgesehen, 
so steigt die Warmeiibergangszahl an der Kondensationsseite bei geeignetem 
Profil der Rillen auf ein Mehrfaches gegeniiber den Werten bei glatten Rohren. 
Der Grund liegt darin, dass die Oberflachenspannungen gewaltige Druckgra- 
dienten in der Ebene senkrecht zur Rohrachse erzeugen. In Figur 2 ist ein 
Rillenprofil dargestellt. In Punkt J ist infolge der Konvexitat der Konden- 
satoroberfliche ein Uberdruck in der Kondensathaut gegeniiber dem Dampf- 
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drucke vorhanden, welcher durch die Oberflachenspannung und den Kriim- 
mungsradius der Kondensathaut @; gegeben ist. Umgekehrt ist im Punkt JJ, 
des Vorzeichenwechsels der Oberflachenkriimmung wegen, ein Unterdruck vor- 
handen. Auf diese Art entsteht ein Druckgradient in der Kondensathaut in 
Richtung von Berg zu Tal. Ware keine Gravitation vorhanden, so wiirde sich 
die Oberfliche der Kondensathaut nach dem umhiillenden Kreis g einstellen. 
Der Unterdruck im Punkt JJ ist bei allen in der Praxis vorkommenden Fallen 
sehr klein gegeniiber demjenigen, welcher durch die «Saughéhe» normal zur 


Fig. 1 Fig. 2 
Vertikales, aussen langsgerilltes Rohr. Kondensatprofil an einer feingewellten Oberflache. 


Bildebene (Figur 2) hervorgerufen wird. Die Oberflache der Kondensathaut 
wird somit in der Talgegend hineingesogen, und zwar in dem Masse, als das 
durch die Gravitation im Tal abfliessende Kondensat die Héhenenergie in Rei- 
bungsenergie umwandeln kann. 

Am Rillenberg ist der Druck des Kondensats infolge der Oberflachenspan- 
nung grdsser als der Dampfdruck. Es muss somit eine Arbeit geleistet werden, 
wenn das Kondensat von der Grenze des Dampfraumes in die Kondensat- 
schicht eindringt. Um eine Idee iiber die Verhaltnisse zu erhalten, sei hier 
trockener, gesattigter Wasserdampf von 100°C und ein Kriimmungsradius am 
Scheitel des Rillenberges von 0,6 mm vorausgesetzt. Damit ergibt sich ein 
Temperatursprung zur Aufbringung der obenerwahnten Arbeit nach den Be- 
ziehungen von Lord KELVIN und von HELMHOLTZz tiber das thermodynamische 
Gleichgewicht zwischen dem Wasserdampf und der gegeniiberstehenden Was- 
seroberflache von nur 1,7 + 10-5 °C. Dieser kleine Temperatursprung wird in 
den nachfolgenden Rechnungen vernachlassigt. 

Um den Vorgang der Hautkondensation an feingewellten Oberflachen besser 
behandeln zu k6nnen, seien die durch weiter unten entwickelte Rechnungen 
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_ bestimmten Gréssen (Hautdicke h, Druck #, Kondensatgeschwindigkeit w und 


Warmeiibergangszahl a) an einem konkreten Beispiel erértert (Figur 3). 
Wie man sieht, erzwingt der grosse Druckgradient eine ausserordentlich 
diinne Kondensathaut (einige um im Rillenscheitel). Daraus geht die Notwendig- 


_ Warmeibergangszah! ce 
ia keal/: Un) 


g/ 42 O3 mm 
=—~ Abgewickelter Unfang des Rillenprofils 


Fig. 3 
Hautdicke h, Druck , mittlere Geschwindigkeit w und Warmeiibergangszahl. 


keit einer Priifung des Einflusses der Wandrauhigkeit hervor. Er wird im Ab- 
schnitt 3.4 besprochen. 

Die Warmeiibergangszahl ist, da es sich um ein reines Warmeleitproblem 
handelt (Reynoldssche Zahlen von der Gréssenordnung 0 bis 10), umgekehrt 
proportional der Dicke der Kondensathaut. Da letztere so klein ist, wird die 
Warmeiibergangszahl sehr gross und erreicht im Beispiel von Figur 3 mit 
100000 bis 150000 kcal/°C m2h Werte, die sogar teilweise iiber denjenigen 
liegen, welche GNaM bei Tropfenkondensation gemessen hat. Dabei ist jedoch 
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zu bemerken, dass bei gerillten Rohren die Kondensation sich nur auf die 
Gegend um den Rillenscheitel beschrinkt, weil das Rillental mit der grossen 
Dicke des Kondensats relativ viel schlechtere Warmeiibergangszahlen auf- 
weist. Die auf die totale Rohroberflache bezogene mittlere Warmetibergangs- 
zahl ist natiirlich viel kleiner als die oben angegebene, jedoch noch immer ein 
grosses Vielfaches derjenigen bei glatten Rohren. | 


3. Berechnung des Warmedurchganges 


3.1. Rechnerische Bestimmung des Profils der Wellen 


Ausser den schon unter 2 angefiihrten Gréssen definieren wir noch: 
h = Dicke der Kondensathaut, normal zur Wandoberflache gemessen, 
o = Kriimmungsradius der Oberflache der Kondensathaut, 
k = Kapillarkonstante (Kondensat gegen Dampf), 
= Druck im Kondensat, konstant tiber 4 angenommen, 
mittlere Geschwindigkeit des Kondensats, 
Kondensatmenge pro Langeneinheit (parallel zur Rohrachse) und 
abgewickelte Lange des Wellenprofils vom Scheitel aus gemessen. 
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Fig. 4. 


Der ausserordentlich kleinen Reynoldsschen Zahlen der Kondensatstré- 
mung wegen wird laminare Strémung vorausgesetzt. Aus dem gleichen Grunde 
wird die Anderung der kinetischen Energie gegeniiber derjenigen der Druck- 
energie als unterdriickbar angesehen. Der Druckverlust ergibt sich somit auf 
der infinitesimalen Strecke ds zu (Figur 4) 


dp = —3 = ds. (2) 
Der durch die Oberflachenspannung gegebene Druck ist der Kriimmung 
der Oberflache proportional. Man kann somit 


k 


rake (3) 
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‘setzen. Der Zuwachs an Kondensatmenge auf der infinitesimalen Strecke ds 
ist durch die Warmeleitungsgleichung 


gegeben. Weiterhin besteht noch die Gleichung 
q=ywh. (5) 


Ferner muss die Form des Wellenprofils gegeben sein. Die Vorschrift dieser 
Form fihrt zu direkt unlésbaren Differentialgleichungen. Sie miissen mit dem 
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Fig. 5 
Profil der den Differenzenrechnungen zugrunde gelegten feingewellten Oberflache 
und das der gefrasten Versuchsrohre (Masse in Millimeter). 
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_ Differenzenverfahren geldést werden. Um jedoch zu einer moglichst giinstigen 
“Form zu gelangen, schlagen wir eine Vorschrift tiber den Verlauf der Konden- 
satdicke h in Funktion des Weges s vor, und zwar in der Form 


h= Ps). (6) 


Gleichungen (2) bis (6) ergeben nach elementaren Operationen: 
1 1 SyhO | / ds | ds 7) 
es) ryk J {2 Dis) | P(g)’ 


-wobei gy, den Kriimmungsradius der Kondensatoberflache fiir s = 0 bedeutet. 
“Wenn man von der Warmeleitung in der gewellten Wand absieht, ist jedes 
_Wegintervall As fiir den Warmeiibergang gleichwertig, mit anderen Worten: 
‘Keine Stelle s; wird vor der Stelle s, bevorzugt. Diese Uberlegung fihrt zu 


Se eee Te ee ee eee 


—wa na eT,e S 
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einem Wellenprofil mit konstanter Dicke der Kondensathaut. Wir setzen somit: 


h = @(s) = ho = const. 


Gleichung (7) wird zu 
1 1 SAS 
Qo ark Ihe 


und erméglicht die Berechnung des Verlaufes der Oberflache des Kondensats 
mittels einer graphischen Integration. Das Wellenprofil ist hiermit auch ge- 
geben, da ja h = hy = const ist. 

Auf die Stoffwerte wird noch in Abschnitt 3.3 eingegangen. 

Urspriinglich wurde beabsichtigt, die gerillten Rohre bei den Kondensatoren 
der Dampfturbinenanlagen zu verwenden. Das durch graphische Integration 
erhaltene, in Figur 5 dargestellte Profil wurde fiir eine Kondensationstempera- 
tur des Wasserdampfes von 30°C und die Werte h) = 0,005 mm, # = 0,53°C 
und @9 = 0,6 mm nach Gleichung (8) ermittelt. Der Bestimmung der Arbeits- 
weise dieses Profils bei anderen Zustaénden diente, wie unten ausgefiihrt wird, 
eine Differenzenrechnung. 


3.2. Berechnung des Warmeiiberganges 


Werden die verschiedenen Gréssen eines bekannten Zustandes fiir s = s,, 
mit dem Index m und an der Stelle s + As mit +1 bezeichnet, so ergeben 
sich gemass den Gleichungen (2) bis (5) und Figur 6 


1 
Ont ™ 1 3 Gnt1 ‘As °) 
a ha Ry 
Ad, = (10) 
Inti = In = eye , (11) 
As ; 
As 
Ay = — 
y ae (13) 
1 
huss = hn + [(P— vy) + Z (Ap — 4y)] As, (14) 


wobei die Bedeutung des Radius R, der Winkel m und wp der Figur 6 zu ent- 
nehmen ist. Das der Differenzenrechnung zugrunde zu legende Wertepaar 
(49; Qo) ist so zu wahlen, dass im ebenen Schnitt senkrecht zur Rohrachse zwei 
einander zugeordnete «tiefste» Talpunkte der als stetig vorausgesetzten Rillen 
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und Kondensatbegrenzung auf ein und denselben Rohrradius zu liegen kom- 
-‘Inen, womit auch der Forderung nach Symmetrie gentigt wird. Fir jeden 
Dampfzustand und Wert der Temperaturdifferenz % kénnen somit verschie- 
dene Profile des Kondensats ermittelt werden. Man erhalt die in Figur 7 dar- 
es Kurvenschar, wobei jede einzelne der Kurven die Gleichungen (9) bis 
q (14) und die soeben besprochenen Grenzbedingungen erfiillt. Jedem von diesen 
_ Profilen entspricht eine gewisse Kondensatmenge der Talstr6mung, die sich 


Fig. 7 


Kondensatprofile bei verschie- 
Fig. 6 denem Fiillgrad der Ablaufrinne. 


Ee NG wee a ee ey ee Te oe  ? 


damit in Funktion der Rohrlainge bestimmen lasst. Die Kondensatmenge, 
welche in der Richtung der Rohrachse bei gegebenem Profil in der Rille ab- 
fliessen kann, wurde durch ein Relaxationsverfahren nach SOUTHWELL be- 
-stimmt. Figur 8 legt das Endergebnis aller Berechnungen der mittleren (auf 
| die dem Radius des Rillentales entsprechendeZylinderflache bezogenen) Warme- 
_tibergangszahlen in Funktion von Rohrlange, Dampfzustand und Temperatur- 
 differenz zwischen Sattdampf und Wand dar. 

| Der plétzliche starke Abfall der Warmeiibergangszahl bei einer gewissen 
_Rohrlange ist durch den Umstand gegeben, dass die Rillen im unteren Teil des 
Rohres vollaufen. Hohe, durch die Kapillarkrafte hervorgerufene Druckgra- 
-dienten werden in diesen Partien vollstandig ausgeschaltet. Fiir die tiber- 
-schwemmten Rillen in den unteren Partien der Rohre ist, wie schon oben er- 
_wahnt, laminare Str6mung angenommen worden, und die Rechnung des Warme- 
iiberganges erfolgte nach NussELtT. Diese Annahme musste getroffen werden, 
da die kritische Reynoldssche Zahl fiir eine Kanalstr6mung mit Profil nach 
Figur 5 bis jetzt noch unbekannt ist. Vermutlich ist, besonders bei heissem 
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Kondensat, in der Rille Turbulenz zu erwarten, womit die wirklichen Warme- 
iibergangszahlen im Gebiete des starken Abfallens (Figur 8) hoher als die hier 
berechneten liegen diirften. 
Andererseits ist es auffallend, wie die Warmetibergangszahl dort, wo die 
Rillen noch nicht vollgelaufen sind, eine sehr schwache Funktion der Rohrlange 
darstellt. Dies ist damit zu erklaren, dass die 6rtlichen Warmeiibergangszahlen 
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Kurve 4 b by bs c 
t, 30 100 100 100 200°C 
o 0,7 10 oy V0) Sc 
Fig. 8 
Mittlere Warmetibergangszahl in Funktion der Rohrlange fiir verschiedene Sattdampftemperaturen 
und Temperaturdifferenzen #. 


in den Scheiteln um so grésser werden, je mehr die wirksame Oberflache durch 
das Anschwellen des Kondensatstromes in der Rille vermindert wird (Figur 7). 
Die Kompensation beider Effekte ist so weitgehend, dass sich mit zanehmender 
Rohrlange nur wenig abnehmende, fast konstante Warmetibergangszahlen er- 
geben. 


3.3. Verdnderlichkett der Stoffwerte 


Von allen Stoffwerten ist die Zahigkeit des Kondensats diejenige, deren 
Veranderlichkeit iiber die Dicke des Kondensatfilms den starksten Einfluss 
auf den Warmetibergang hat. Eine Berechnung der mit der Filmdicke linear 
verdnderlichen Zahigkeit ergab jedoch, dass man fiir in der Technik tibliche 
Betriebszustande einfach mit dem arithmetischen Mittelwert der Zahigkeiten 
an der Dampf- und Wandseite rechnen kann. 
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3.4. Einfluss der Oberflachenbeschaffenheit 


. Die Berechnung des Warmeiiberganges bei einer sinusférmig angenomme- 
nen Wandrauhigkeit zeigt, dass bei Rauhigkeiten von der Gréssenordnung 
_< 0,2 Filmdicken der Einfluss derselben vernachlassigbar ist. 


3.5. Warmefluss durch die Wand 


a 
ll 
4 Die Berechnung des Warmeflusses durch die Wand (Methode der Quadrate) 
_ fihrte zur Erkenntnis, dass man in Fallen von warmedurchgangsmissig hoch 
4 geziichteten Kondensatoren die Wand des Warmedurchgangs aus gut leiten- 
4 dem Material herstellen muss. Gelten die Bezeichnungen: 
j Ay» = Warmeleitfahigkeit der Wand, 

_ da = minimale Wanddicke (1 mm in bezug auf die Abmessungen von Figur 5), 
4 
4 


Q = iibertragene Warmemenge, 
F 


I 


Oberflache des dem Rohr koaxialen Zylinders entsprechend dem Rillen- 
tal und 

Ai,,= mittlerer Temperatursprung in der Wand, 

so ergibt sich fiir alle der Figur 5 geometrisch ahnlichen Wandquerschnitte 
254. 20 


Aip= (15) 


4. Messungen 


3 Um die theoretischen Ausfiihrungen tiberpriifen zu kénnen, miissen natiir- 
. lich die der Rechnung zugrunde gelegten Voraussetzungen moglichst weit- 
_ gehend erfiillt werden. Vor allem gilt dies fiir Grésse und Form des Profils der 
_ Rillen. Die kupfernen Versuchsrohre von 525mm Lange, 16 mm Innen- und 
- 19,2 mm minimalem Aussendurchmesser wurden durch Frasen am Umfang mit 
- 29 Rillen versehen. Die Toleranz der Profilmasse betrug + 0,01 mm. Das 
_ Rillenprofil ist in Figur 5 dargestellt. Bei den Versuchen wurde Sattdampf von 
- 100°C verwendet. Die Bestimmung der Warmemenge erfolgte durch Konden- 
satmessung. 

4 Versuche in den Escher-Wyss-Laboratorien, Ziirich, wurden durchgefithrt : 
| a) an glatt polierten Rohren, welche kurz vor dem Versuch 48 h in Létwasser 
- __ gelegt worden sind (um Hautkondensation zu erzwingen), 

_b) an ebenfalls glatt polierten, jedoch noch (zur Erzielung guter Benetzung) 
mit dem Waschmittel «Vel» behandelten Rohren und 

c) an oxydierten Rohren ohne jegliche Behandlung. 


: Durch Messung des Kondensats, der Temperatur des Sattdampfes, des 


_ Kiihlwassers und der Oberfliche der Versuchsrohre wurde die Warmedurch- 
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gangszahl bestimmt. Um die Warmeiibergangszahl auf der Dampfseite be- 
rechnen zu kénnen, musste, ausser Gleichung (15), noch die Warmetibergangs- 
zahl auf der Wasserseite bekannt sein. Diese wurde nach der Beziehung von 
KRAUSSOLD bestimmt. Um den Einfluss der Ungenauigkeit dieser Formel mog- 
lichst klein zu halten, wurden in den Rohren sehr grosse Wassergeschwindig- 
keiten eingestellt. 
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10000 * 
2 J OST Woe Fe eh EM a; COPECO, 
a = gerechnete Kurve fiir Rauhigkeit = 0, O ohne Oberflachenbehandlung, 
* Oberflache mit L6twasser behandelt, A Oberflache mit «Vel» behandelt. 
Fig. 9 


Versuchsergebnisse an gefrasten Versuchsrohren. 


In Figur 9 sind die Ergebnisse von Versuch und Rechnung dargestellt, 
woraus man die gute Ubereinstimmung zwischen den beiden Resultaten er- 
sehen kann. 

Gefraste Rillenrohre kommen natiirlich fiir die praktische Anwendung nicht 
in Frage, da sie zu teuer sind. Nach dem guten Ergebnis der Versuche an genau 
gefrasten Rohren dachte man vielmehr an gezogene Rohre. Bei diesen kénnen 
aus Herstellungsgriinden die Masstoleranzen nicht so klein gehalten werden 
wie bei den ersteren. Die Versuche wurden deshalb an einem 2510 mm langen, 
gezogenen, gerillten Rohr wiederholt, um den Einfluss der Ungenauigkeit des 
Rillenprofils auf die Warmeiibergangszahl zu tiberpriifen (Figur 10). 

In Figur 11 ist die Warmetibergangszahl an der Kondensationsseite in Funk- 
tion der Temperaturdifferenz zwischen Dampf und Wandtemperatur aufge- 
tragen. Die Streuung ist auf die bei unbehandelter Oberflache durch ein Schau- 
glas beobachtete partielle Tropfenkondensation zuriickzufithren. Die gestri- 
chelte Kurve stellt die Rechnungsergebnisse dar; die ausgezogene hingegen ist 
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Fig. 10 


° Sappeler-Chromséure uComprax 


6 Chromséuré 


Fig. 11 
Versuchsergebnisse an gezogenen Versuchsrohren. 


: die untere Umhiillende der Messergebnisse. Das Resultat erscheint im Hinblick 
auf die praktischen Anwendungen als vollauf befriedigend. 


5. Modellahnlichkeit 


Um sich die Arbeit der Versuche mit Dampfen anderer Fliissigkeiten zu 
ersparen, werden hier einige aus einer Ahnlichkeitsbetrachtung hervorgegan- 
_ gene notwendige und hinreichende dimensionslose Kenngréssen angegeben. 
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5.1. Laminare Strémung in der Ableitungsrille des Kondensats 


Die dimensionslose Schreibweise von Gleichung (9) und der Gleichung des 
Druckverlustes in der Ableitungsrille fiihrt auf folgende zwei dimensionslose 


Kenngrossen 


y. 
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Graphische Auftragung der Warmetibergangszahl mittels dimensionsloser Kenngréssen. 


Bei laminarer Str6mung des Kondensats in den Ableitungsrillen wird die 
Nusselt-Zahl somit eindeutig durch 


ee eae Wt os and 
Nu = A eee at ea) 


gegeben (siehe Figur 12). 


5.2. Turbulente Stromung in der Ableitungsrille des Kondensats 


Mit dem Ansatz von Brasius fiir den Widerstandsbeiwert und diesmal auch 
bei Beriicksichtigung des Warmeitiberganges in der Ablaufrille hat man 


214 QT/4 lla 
Nu = po( HAD ~~ qild OIL 


TRY Oy’ O84 p27 gale ? 


Re; Pr) (16) 


aoe 2S. YG 


=. Sa Tes 4” 
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mit - 


a B Oy 
yng 


und Pr = ene 


Anstatt der zwei Variablen (Re; Pv) kann man erfahrungsgemass 
Rew py? 


mit etwa m = 3/4 und = 1/3 schreiben, womit sich die Anzahl der Variablen 
in (16) von vier auf drei reduziert. 


6. Zusammenfassung 


Nach einer einleitenden Ubersicht der Haut- und Tropfenkondensation wird 
die Hautkondensation an einer unebenen Flache erértert. Die Oberflachen- 
spannungen rufen, der verschiedenen Kriimmungen der Kondensatoberflache 
wegen, sehr starke Druckgradienten hervor. Die Wirkung letzterer kann ein 
Vielfaches derjenigen des Schwerefeldes betragen (bis zu einigen Zehner- 
potenzen). Starke Druckgradienten ergeben aus Gleichgewichtsgriinden diinne 
Kondensathaute, was nach der Theorie von NUSSELT zu sehr grossen Warme- 
tibergangszahlen fiihrt. Fiir einen Kondensator mit achsenparallel gerillten, 
vertikalen Rohren werden bei gegebenen geometrischen Verhaltnissen die durch 
Differenzenrechnung erhaltenen Warmetibergangszahlen angegeben. Letztere 
wurden durch Versuche bestatigt. Kennzahlen der Modellahnlichkeit bestimmen 
die Nusseltsche Zahl des Warmeiiberganges. 


Summary 


Following a synopsis on film- and drop-condensation, one deals with con- 
densation on an uneven surface. The surface tensions give rise to very great 
pressure gradients because of the different curvatures of the surface of the 
condensate. The action of the latter can amount to a multiple (up to a few powers 
with basis ten) of the value corresponding to the gravitation. Large gradients of 
pressure produce thin films of condensate, and according to the theory of NUSSELT 
they lead to large coefficients of heat transmission. Heat transfer coefficients are 
given determined with calculus of finite differences for a condenser with vertical 
tubes, grooves parallel to the tube axis and given geometrical form. The latter 
are confirmed by tests. Adimensional numbers determinate the NussELT number. 


(Eingegangen: 19. Februar 1953.) 
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Die exakte Losung der Integralgleichungen 
gewisset Schwingungsprobleme 


Von Jens RarneR Marta Rapok!) und ALFRED HELLER, Zirich*) 


1. Einleitung 


Die Schwingungsprobleme gewisser dynamischer Systeme mit polynomialen 
Massen- und Steifheitsverteilungen lassen sich leicht auf Fredholmsche Integral- 
gleichungen mit entarteten symmetrischen Kernen zuriickfiihren. Die Lésung 
solcher Integralgleichungen hat jedoch bisher grosse Schwierigkeiten bereitet, 
so dass man fast immer gezwungen war, Naherungsmethoden zu verwenden. 

In dieser Abhandlung ist eine Methode entwickelt worden, die mit dem be- 
kannten Iterationsverfahren®) verwandt ist. Es wird namlich gezeigt, dass man 
durch geschickte Wahl der Ausgangsfunktion, die durch den Kern bestimmt 
wird, zu Potenzreihen fiir die Eigenfunktionen kommen kann. Die Koeffizienten 
dieser Potenzreihen geniigen besonders einfachen Rekursionsformeln, mit deren 
Hilfe man schliesslich nur eine endliche Anzahl willkiirlicher Konstanten erhalt. 
Diese restlichen Konstanten fiihren dann auf die charakteristische Gleichung 
fiir die Eigenwerte. 

Es liegt in der Natur des eben angedeuteten Prozesses, dass der Ubergang 
zu den exakten Lésungen und der exakten charakteristischen Gleichung rein 
formell durchgefiihrt werden kann. Die Potenzreihen, die man durch diese 
Methode erhalt, haben Koeffizienten, die Potenzen der Frequenzen im Zahler 
und Fakultaéten im Nenner enthalten, so dass in jedem Fall die Konvergenz- 
betrachtungen sehr einfach werden. 

Im zweiten Teil dieser Mitteilung wird die Integralgleichung fiir symme- 
trische Biegeschwingungen von Balken abgeleitet, wobei die Kriimmung der 
Balkenachse als unbekannte Verdnderliche eingefiihrt wird. Im Falle dieses 
einfachen Problems bereitet die direkte Ableitung dieser Gleichung aus der 
Differentialgleichung keine Schwierigkeit. Die asymmetrischen Biegeschwin- 
gungen lassen sich genau gleich behandeln, wahrend es in komplizierteren Fal- 
len vorzuziehen ist, zum Beispiel das Hamiltonsche Prinzip anzuwenden’). 


1) Sonst Aeronautical Research Laboratory, Department of Supply, Melbourne, Australia. 


) 
2) Institut fiir Flugzeugstatik, ETH. 

3) W. ScHNEIDLER, I ntegralgleichungen mit Anwendungen in Physik und Technik (Akademische 
Verlagsgesellschaft, Leipzig 1950). 


4) J. R. M. Rapox, Dynamic Aeroelasticity, Rep. 58 (College of Aeronautics, Cranfield, Eng- 
land). 
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Die allgemein zu verwendende Methode ist im dritten Teil fiir einen spe- 
ziellen Kerntypus entwickelt, wahrend in den letzten zwei Teilen diese Methode 
auf die Lésung zweier Probleme angewandt ist, deren exakte Lésungen seit 
langem bekannt sind. Das erste dieser Probleme befasst sich mit den Schwin- 
gungen eines gleichférmigen Balkens!), wahrend das zweite diejenigen des fest 
eingespannten Keiles?) untersucht. In beiden Fallen wird gezeigt, dass die 
Methode auf die exakten Lésungen fiihrt. 


2. Ableitung der Bewegungs¢gleichungen 


Es ist bereits in der Einleitung erwahnt worden, dass die Untersuchung 
dieser Abhandlung auf der Basis der zustandigen Integralgleichungen ausge- 
fiihrt werden soll. Wie sich sogleich zeigen wird, sind diese Gleichungen von 
demselben Typus, unabhangig davon, ob der Balken an einem Ende einge- 
spannt oder ob beide Enden frei sind. Nachdem die Integralgleichung in jedem 
Falle die Randbedingungen implizit enthalt, das heisst, nachdem ihre Loésungen 
automatisch diese Bedingungen erfiillen, wird es das Ziel dieses Paragraphen 
sein, eine allgemeingiiltige Integralgleichung fiir die zu betrachtenden Falle 
abzuleiten. ; 

Man betrachte den Balken von Figur 1, der am Ende y = 0 eingespannt 


Figur 1 


und am Ende y = / frei sei. Die Massenverteilung und Steifheit des Balkens an 


der Stelle y seien 
m(y) und EI(y) , 


so dass bei dynamischer Belastung, auf Grund der Bernoulli-Euler-Formel, die 
Gleichgewichtsbedingung an der Stelle y folgende Form hat: 


m(y) 22 4 25 (ET) Sa) — my) + (E10) 27" =0. (2.1) 


1) S. TimosHENKO, Schwingungsprobleme der Technik (Springer, Berlin 1932). 
2) G. R. Kircuxorr, Uber die Transversalschwingungen eines Stabes von veradnderlichem Quer- 


schnitt, Berliner Mber., 1879. 
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Dazu kommen die Randbedingungen 


2(0,t) =.2'(0, t) = E I(y) 2"(y, ),-. =[E L(y) 2", O],-1=9, (2. 2) 


: 


von denen die zwei letzten das Verschwinden des Biegungsmomentes und der © 


Querkraft am freien Balkenende ausdriicken. 

Integriert man die Gleichung (2. 1) zweimal von u bis / und von v bis J, so 
erhalt man, unter Beniitzung der zwei letzten Bedingungen von (2. 2), die 
Gleichung 


de aie 


1 L 
02 
[au f my) Se dy + E It) eae ig (2. 3) 


die sich jetzt auf das Gleichgewicht der Biegemomente im Abschnitt y = v und 
y =I bezieht. Die Gleichung (2. 3) ist aber noch keine reine Integralgleichung, 
indem ihre Lésungen jetzt zwar die letzten beiden Bedingungen von (2. 2) er- 
fiillen, die bei ihrer Ableitung beniitzt wurden, nicht aber die tibrigen Bedin- 
gungen (2. 2). 

Um diese Bedingungen auch noch in Betracht zu ziehen, kann man jetzt 
auf zwei verschiedene Arten verfahren, von denen die eine auf die Integral- 
gleichung fiir die Auslenkungen z, die andere auf eine 4quivalente Gleichung 
fiir die Kriimmungen 2” fiihrt. Die zweite Methode soll hier verwendet werden, 
da dann die Lésungen direkt zu Biegemomenten fithren, die vom praktischen 
Standpunkt aus von grésstem Interesse sind. Dazu kommt, dass man im Falle 
einer Naherungslésung die Auslenkungen, das heisst die Schwingungsformen, 
durch zweimalige Integration erhalten kann, wodurch etwaige Fehler bekann- 
terweise reduziert werden. 

Zur weiteren Transformation der Gleichung (2. 3) betrachte man jetzt die 
folgende Beziehung, die auf Grund der ersten zwei Bedingungen von (2. 2) be- 
rechtigt ist: 


. y 0% we Ad) 0 ay oe 7 
Hy,t) = | 5 dn = fan [Ser at = fan [2G Has. 24) 
0 0 0 0 0 


Bevor man von (2. 4) in (2. 3) einsetzt, sollen die Doppelintegrale in diesen 
beiden Ausdriicken umgedreht werden. Unter Bentitzung der Figuren 2a, 2b 
findet man leicht an Stelle von (2. 3) und (2. 4) 


1 


| m0) 20.9 (y — 2) dy + ET) 2"(0, 1) =0 (2. 3") 
und : 


Vv 


fPCNY-H a =%0,0. 2.4’ 


0 
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Setzt man jetzt (2. 4’) in (2. 3’) ein und reduziert wiederum das Doppelintegral 
durch Vertauschen der Integrationen auf ein einfaches Integral (unter Beniit- 
zung von Figur 2c), so erhalt man schliesslich die endgiiltige Integralgleichung 


ila la 


Figur 2a met 2b Figur 2c 


RTE CE eye ny Ee ee a ee 


yr en eo 


_ fiir die Schwingungen des eingespannten Balkens, wobei die Kriimmung 2” die 
unbekannte Veranderliche ist: 
L 
[HG #69 do + BT) 20,1 =0 (2. 5) 
0 
mit 


. sa) = [my (u— Cl) (u— y) du, (2. 6) 


Dabei steht oy fiir den grésseren Wert von € oder y. 

Als nachstes betrachte man den sich frei im Raum bewegenden Balken. Wie 
bereits in der Einleitung erwahnt wurde, soll die augenblickliche Untersuchung 
auf symmetrische Verformungen des Balkens beschrankt werden. Beim freien 
Balken verlangen Impuls- und Drallsatz, dass der Gesamtimpuls und der Ge- 
samtdrall dauernd verschwinden. Aus dem Impulssatz folgt zunachst, dass der 
Schwerpunkt S des Balkens, dessen Lange nun 2/ sein soll, sich nicht bewegt. 


oe 


Z 
Figur 3 


In dem mitbewegten Koordinatensystem (y, z) (siehe Figur 3) besitze S die 
Koordinaten (0, z)). Die Auslenkung des Balkenelementes an der Stelle y be- 
ziiglich des ruhenden Schwerpunktes ist dann 


2(y, t) — 2o(y, #) 
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und seine absolute Beschleunigung 
aly, t) = 2o(s t) = 
Die Bedingung, dass die zeitliche Anderung des Gesamtimpulses verschwinde, 
nimmt daher folgende Form an: 
+1 
[my 20,9) dy — 2M 4(t) = 0, (2.7) 
= 


wobei 
M = Masse des halben Balkens. 


(2. 7) kann man jetzt wieder so umformen, dass unter dem Integral z durch 4 
ersetzt wird. Die Gleichung (2. 7) erhalt dann die Form: 


+1 


wX(C) 2"(C,t) dl — M2,(t) =0, (2.7') 


Wo 
= [ mo) (v — C) dv. (2. 8) 
(s 


(Der Drallsatz liefert bei den symmetrischen Schwingungen keine Aussage.) 

Um jetzt die Integralgleichung fiir den freien Balken abzuleiten, kann man 
gleich vorgehen wie beim einseitig saben Stab. Die Randbedingungen 
(2. 2) behalten ihre Giiltigkeit. In (2.1), (2.3), (2.:3’) ist an Stelle von z die 
neue Beschleunigung z — 2) einzusetzen. Statt (2. 5) erhalt man so: 


[700393 Ct) al — wX(y) Zl) +E T(y) 2"(y,)=0. 2.5") 


Eliminiert man jetzt z) zwischen den Gleichungen (2. 5’) und (2. 7’), was im 
Sinne der allgemeinen dynamischen Theorie bedeutet, dass z) eine verborgene 


Koordinate ist, so findet man schliesslich die gesuchte Integrodifferentialglei- 
chung 


l 
E Ty) 2"(y,t) +] Kly, 2", 8) ae =0, (2. 9) 
0 
wo 
> Ss = a L(y 1 
Kiy cy yee ee (2. 10) 


Bisher wurde angenommen, dass die Masse des Balkens stetig verteilt sei. Man 
sieht aber aus der Definition der Gréssen uw, 2, dass man in der Balkenmitte 
eine konzentrierte Masse 2 My hinzufiigen kann, ohne diese Massenfunktionen 


mee 


a 


7 


tee CAS Oe 2 eh NUNS ee ee a 
bias “ at . Aa 


i a ae he 


wey) oe 
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zu andern, da sie nur von der Masse zwischen der Schnittstelle und dem freien 
Balkenende beeinflusst werden. Definiert man also die Gesamtmasse des halben 
Balkens mit beliebiger longitudinaler Massenverteilung m/(y) als 


1 
M = | miy) dy+M,, (2. 11) 
0 


so sieht man sofort, dass die Gleichung (2.10), wie zu erwarten war, fiir alle vor- 
geschlagenen Probleme zustandig ist. So erhalt man zum Beispiel fiir My = co 
die urspriingliche Gleichung (2.5) fiir den eingespannten Balken, wahrend 
M, = 0 dem Balken ohne konzentrierte Masse entspricht. 

Schliesslich soll noch eine kurze Bemerkung iiber die Ableitung von ahn- 
lichen Integralgleichungen fiir kompliziertere Probleme derselben Art gemacht 
werden. So mag es zum Beispiel wiinschenswert sein, nicht nur symmetrische, 
sondern auch antisymmetrische Schwingungen zu studieren. Oder es mégen 
Torsionsschwingungen von Interesse sein. Beim Flugzeug handelt es sich meist 
nicht nur um einen einfachen Balken, sondern um zwei oder mehrere Balken, 
die aufeinander senkrecht stehen (Rumpf, Fliigel, Schwanzflache) und Biege- 
wie auch Torsionsschwingungen unterliegen. Alle diese Probleme bedeuten eine 
Verallgemeinerung des oben behandelten Problems in dem Sinne, dass die ela- 
stischen Massensysteme eine gréssere Anzahl von Freiheitsgraden haben oder 
dass noch weitere verborgene Koordinaten vorhanden sind. Den Fall antisym- 
metrischer Schwingungen des freien Balkens kann man zum Beispiel durch 
die Einfiihrung einer Rotation um den Schwerpunkt studieren. Zu den Bedin- 
gungen (2. 2) hat man jetzt aber die Bedingung des Verschwindens des Dralls 
hinzuzufiigen. 

Im allgemeinen ist es aber in komplizierteren Fallen nicht mehr so einfach, 
die zustandigen Bewegungsgleichungen durch ahnliche Uberlegungen abzu- 
leiten. Es ist daher besser, auf die Prinzipien der allgemeinen Dynamik, das 
heisst auf die Lagrangeschen Gleichungen oder auf das Hamiltonsche Prinzip 
zuriickzugreifen. In der Tat sind die Gleichungen auf diesem Wege fiir ein Flug- 
zeug mit Pfeilfliigel und Rumpf in der Arbeit von RApDOK!) abgeleitet worden, 
und es ist dort gezeigt worden, wie die Probleme dieser Abhandlung als Sonder- 


fille herauskommen. 


3. Probleme natiirlicher Schwingungen und ihre angenaherte Loésung 


Setzt man jetzt voraus, dass der Balken natiirliche Schwingungen ausfihrt, 


so kann man schreiben 
z"(y, t) = ert Z (9). (381) 


1) J.R.M. Ravox, Dynamic Aeroelasticity, Rep. 58 (College of Aeronautics, Cranfield, England). 
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Damit erhalt (2. 9) die Form 
U 


EI(y) Z(y) — 08 | K(y, 9) Zp) 4p = 0. (3.2) 


0 


Diese Gleichung wird leicht in eine Fredholm-Integralgleichung zweiter Art 
umgewandelt, indem man folgende Transformation einfiihrt: 


EE TOl Cink ee ee ae 


VE I(y) VE I(9) 


so dass man jetzt anstatt (3. 2) die Gleichung 


Z19)— 08 [ Roo 2 9) dp = 0 (3. 2’) 


erhalt. 

Dieses ist die Grundgleichung der spateren Untersuchungen, und die anzu- 
wendende angenadherte Lésungsmethode soll jetzt kurz beschrieben werden. 
Es soll angenommen werden, dass der Kern K ausgeartet ist in dem Sinne, 
dass er ein homogenes Polynom n-ten Grades in y und y ist). 


PrlQ, Y) = An Y" + Ana V1 PH $a VG" + ay GO" (P<), | aan 
EY, Q) = On Pr Ay GY + 2 ay V1 FE ay VY? ONG) | 


Es ist leicht einzusehen, dass K immer symmetrisch ist, so dass 
LFA: Whesg = ln Or Oles. 


Ein solcher Kern entspricht einer Massenverteilung, die selbst durch ein Poly- 
nom dargestellt werden kann. 

Weiterhin soll die Gleichung (3. 2’) noch normalisiert werden, das heisst, 
das Integrationsintervall soll von 0 bis 1 gehen, so dass die Gleichung jetzt die 
folgende Form hat 


Z(y) — K y, p) Z(p) dp =0. Bey 


Ist die Lange / nicht gleich 1, so verursacht die entsprechende Umformung nur 
eine Anderung des Frequenzparameters w2. 


1) In den Beispielen treten etwas allgemeinere Kerne auf, die aber ebenfalls entsprechend der 
folgenden Theorie behandelt werden konnen. 
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Man betrachte jetzt den Integraloperator 


; 

far 

3 j 

: = | Pale. 9) 2) dp + [20 9) dp , (3. 5) 
“a oo 

Z(p) = y™ 

- m — aya Nan (m — 2) ay_ 

eee ferme rary Mee fn 

a 

: fy (2—n) a, (—n) ay 

4 +more + waaay} | (e5 


en a eee 


GZ 
LS omar es mtn mene L 


m+ 2 


Der Operator transformiert also eine beliebige Potenz von y in die Form 


PPT CROPS 


Lom = yimt0*"C, + Pay), (3. 6' 


- wobei das neue Polynom P%(y) wiederum vom Grad n ist. Dieses Resultat 
gibt einen Hinweis fiir den Ansatz der Potenzreihe, mit der Fredholm-Gleichnu- 
_ gen des obigen Typs immer annahernd gelést werden kénnen. 

Man setze namlich die Lésung Z der Gleichung (3. 2”) in der folgenden all- 
_ gemeinen Form an: 

k 


= S Ay; pi (mt) (3. 7) 


j=0 1/=0 


- Bricht man jetzt die Reihe im Sinne einer Approximation bei 7 = k ab, so ist 
der Fehler 


n 
De DS Ay gy PAD tDt8 Com gayer (3. 8) 


i=0 


Durch Koeffizientenvergleich der Potenzen von y ergeben sich Beziehungen 
zwischen den Konstanten A,. In der Tat fiihrt (3. 7) nach Anwendung von 
_ (3. 6’) bei Vergleich der Potenzen von y mit Exponenten > zur Rekursions- 
_formel 


Ay; = 0 Ajj Ciaymsners 9 =1,2,...,2; 2=0,1,2,...,m) (3.9) 


durch welche man die A, ; durch die A;,) ausdriicken kann. Zudem sieht man 
von (3. 6) und der Definition von C,,,, dass die Koeffizienten A;,; durch A; 9 
-ausgedriickt im Nenner Fakultéten von der Ordnung 7(” + 1) + / enthalten 
werden, wahrend man im Zahler Polynome von w? vom Grad m hat. Daher 
werden die angesetzten Lésungen fiir gentigend grosse k gut konvergieren. 
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Schliesslich gibt der Koeffizientenvergleich der Potenzen von y mit Expo- 
nenten <n das folgende System von ” + 1 Gleichungen: 


a 
2 0 
Ano = ot YS Aus i(@rieret 


a 
2 1 
Avsa= ot 33 Ass j(m+i1)+l4+2 


G=01= 
OR Ae are eam vn Go o> Qabem (3. 10) 
k n ws 
= a2 7 n—1 
Ay, = w? J’ Ay; j(n+ij)+lin ”’ 
7=0 1=0 
k n . 
ens j n 
Ago. =O 2, dy Au Fee le met’ J 


das mit Hilfe von (3. 9) so reduziert werden kann, dass es nur die + 1 Kon- 
stanten A;,, enthalt. Nach dieser Reduktion kann man dann aber die Koeffi- 
zienten A,,9 — Ay,9 durch Ao,» ausdriicken und in der letzten Gleichung ein- 
setzen. Auf diese Weise erhalt man ein Polynom in w? von der Ordnung kn. 
Damit die Konstante A»,9 endliche Werte haben kann, muss w? eine Wurzel 
dieses Polynoms sein. Somit ist dieses Polynom die charakteristische Gleichung 
fiir die Bestimmung der Eigenwerte, und die Eigentanikaoney sind durch 
(3. 7) gegeben. 

Es wird nun von Interesse sein, diese Lésungsmethode an praktischen Pro- 
blemen auszuprobieren. Aus der vorausgehenden Betrachtung ergibt sich ném- 
lich, dass die Wahl einer gréssern Summationsgrenze die Form der Lésungen 
nicht andert. Aus diesem Grunde und wegen der guten Konvergenz der Kon- 
stanten Aj; kann man jetzt vermuten, dass die angenaherten Lésungen (3. 7) 
die ersten Glieder der Potenzreihenentwicklungen der exakten Lésungen repra- 
sentieren. Daraus wiirde folgen, dass die ersten » +1 Glieder des Polynoms 
fiir die Bestimmung der Frequenzen den Gliedern der exakten charakteristi- 
schen Gleichung entsprechen werden. 

Somit erhalt man die Moglichkeit, die Eigenfunktionen und Eigenwerte 
numerisch beliebig genau zu approximieren durch eine entsprechende Wahl der 
Summationsgrenze k. Auf der andern Seite sieht man auch, dass, nach Beob- 
achtung der Gesetzmiassigkeit in der Bildung der Rekursionsformel und der 
Koeffizienten des charakteristischen Polynoms, der Schritt von k zu k+1 
rein formell ausgefiihrt werden kann, so dass die Methode bei numerischen 
Rechnungen den gleichen grossen Vorzug hat, den auch die Iterationsmethode 
bietet, idem jeder Schritt zu einer verbesserten Lésung fiihrt. Durch die 
jetzige Methode erhalt man alle gewiinschten Eigenfunktionen und Eigenwerte 
auf einmal, wobei Anzahl und Genauigkeit durch k festgelegt sind. Vergréssert 
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man k, so erweitert sich die Sakulargleichung fiir die Eigenwerte ohne Ande- 
tung der bereits berechneten Glieder um zusitzliche héhere Potenzen. 

In den folgenden zwei Paragraphen soll die eben abgeleitete Methode auf 

wei Schwingungsprobleme angewandt werden, von denen exakte Lésungen 
bekannt sind. Es handelt sich dabei um die Probleme des schwingenden, ein- 
_ gespannten, gleichmassigen Balkens und des eingespannten vollen Keiles. Dabei 
-sollen jedoch nicht die oben abgeleiteten Formeln beniitzt werden. An deren 
Stelle werden die gewiinschten Formeln direkt abgeleitet, und es wird gezeigt 
werden, dass in der Tat die versprochenen Anndherungen resultieren. 


5 3 


4. Schwingungen von Balken konstanter Massen- und Steifheitsverteilung 


Obwohl am Ende nur der Fall des eingespannten Balkens behandelt werden 
soll, werden doch die allgemeineren Formeln fiir den freien Balken mit einer 
Masse in der Mitte abgeleitet. Setzt man also 


RS eR ee yee nn = eee Ee Qos a8 3. 


m(y)=—-, My=kM, (4a) 


4 so erhalt man von (2. 10), (2. 8), (2. 6) 


eek. 25 Dy) 10) 
F 7 Ky, 6) = L(y, SS i (4. 2) 
wo 
Ty, ¢) = [ (w—¢) (w— 9) dw, (4. 3) 
7,6 
. Ay) = | (w—y) 4 (4. 4) 
3 y 
Daher hat man fiir den Kern 
ie pe eS a i a 
Pe a ee | > 0 | 
l opi-o , f-y Se 
@-9)|—4—-+-5"-Fayl] V<o- 
-Normalisiert man jetzt durch Einfiihrung der Koordinaten 
Lee BS 2 (4. 6) 


, l Si) 
und fiihrt die neuen Veranderlichen 
p=1-Y, p=1-X (4. 7) 


ein, was einer Verlagerung des Koordinatenursprunges an die Balkenspitze ent- 
-spricht, so nimmt die Integralgleichung (3. 2’) die folgende Form an: 


agi IY fs BS pera ok 


+fels , Hf cate] 2 4, 


(4. 8) 


60: Jens Ratner MartA Rapoxk und ALFRED HELLER ZAMP 


wobei 
2 3 
ie ee _ (K= EI = Balkensteifheit) (4. 9) 


Den Integraloperator auf der rechten Seite transformiert nun die Potenz y” 
in den Ausdruck 


Gira) rh gobs 1 1 R 
(n+ 4)! 6(m+1) Mpa a ea Sie 


Daher kann man das allgemeine Resultat von m Iterationsschritten mit der An- 
fangsfunktion 
Ay + By, 


die durch die Formel (4. 10) angedeutet ist, in der Form 
= de * . ? 
Zigh es DT {A wits ae Bags tary (4. 11) 
j=l 
schreiben. Postuliert man nun den Fehler 


grit mo Ih ye pada ee oy 


4Am (4 m + 3)! ve (4 m + 2)! : 


(4. 12) 


so erhalt man durch Koeffizientenvergleich der Potenzen von m nach Anwendung 
von (4. 10) 


As-1 31 M-1 2! 
A; = a Sareea 1> (4. 13) 
m 4; : 
Ath aa AD erage | 
BA ee : \ 
AI Saapey sasneiea, | ee 
=] wae. : ds 
a aSy  4dl + kag ay] oe 


wobei die letzten zwei Bedingungen von gy? und g? herriihren. Ersetzt man in 
(4. 14) A,,.B; durch ihre Werte (4. 13), so findet man schliesslich die zwei homo- 
genen Gleichungen fiir A,, B,, die eine endliche Lésung dann und nur dann 
haben kénnen, wenn die folgende Determinantengleichung erfiillt ist. 


m 1 wm ai 
1+ : s 
ail 3 Pee to aye 
ike a { 1 m yi 1 
3 aes Se es © ! ine D. = 0 
AGi-y worHayey 34 a, (47-2)! pEeH TED 
“ota 1 
CEs E are 


(4. 15) 
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3 Jetzt fiihre man die Voraussetzung ein, dass der Balken eingespannt ist. Zu 
diesem Zwecke setzt man k = oo und erhalt aus (4. 15) das Polynom in A: 


m m 


eeAfa ad Xe i yi ote 
1 ai os = 
7 on esses mm (47)! 2: aT Dial | 


baa) 


ai 


j= l=1 


jab 


4. 16) 


m 


, ( 
ye a WU 
a ear EY CPSU is | 


j=1 l= 


xe 


Da aber 


ee te 


! 1 g—45 
(4g 4 9) {49 = 20 (aay 


m qu a ; 7 1 1 
rast Gam (47—1)!47 age varias ayy? (42 — 2)! 49 


2m—1 2m-1 


; 1 
Oe za e [4G —le m4 (—m+1)—2]!4(.-me])’ 


=m 


— ve) = i! 
& GF aTay LiCTSTr 


q t=1 


s 


PE EROS EE YR SN ee Oe eT ee eye 
t ¥ }3 


a 
, 


pa 4 G—t44)—11)Gl—1)l (41 £4) 


a 


il i 


= 
2m-1 2m—-1 
1 


4 BO 2) agate m I= we 


y=m 


—e 


7 nimmt (4. 16) die folgende Form an: 
4 m1 j 1 
zt 2 w| 3 [oy — bor 1) 47 4 at 
| are 


(4G — 74 MAI nei+ i - Wa 


2m—-1 2m—-1 


4 “ 1 
5 jit ; 
7a Da [4G—l4+ m)]! [4 7—m+i)4+ 1! 


{[4(l— m+] [4 (7 -—14+ m)] —[4 (I-m+1)-1)[4(lL-m+1) +1} =0. 
(4. 16’) 


‘Dieser Ausdruck kann aber noch vereinfacht werden. Dieses soll aber nur mit 

dem Glied in eckigen Klammern unter dem ersten Summenzeichen durchgefiihrt 

werden, um zu zeigen, dass die ersten Glieder des Polynoms mit den entspre- 

chenden Gliedern der exakten charakteristischen Gleichung iibereinstimmen. 
-Fiir den eingespannten Balken ist diese bekannterweise 

6 475 ay 2a et (—1)7+1 ;| 

E 4 = 1 ——_.-—___— A }=0. 4.17 

cosh V4 cos VA + 1 = 2 +2, (4G+ i)! | ( ) 
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Somit muss die folgende Identitat bewiesen werden, da man fiir 7 = 0 unmittelbar 
das Glied —A/12 erhalt: 


LAT 4 aed Pe ee 249 ER ee 
_ ‘ eG Te MTD! Tra” eiea 


Nun ist aber 


3 [4 7 Grae {4 Gb Ud ee 
Sie ate Gel -G get od 
=o (es le te ake 


da 
Schliesslich hat man 


Soiled bea Cop a Eee epee ce = 


t=1 


l| 


=47+2+4 8 ate 4) Si+4 + (1— q) tite} =4j71+24 22741 (—1)3+1, 


was zu beweisen war. 

Man sieht also, dass die (m+ 1) ersten Glieder der charakteristischen Glei- 
chung (4. 16’) mit den entsprechenden Gliedern der exakten Gleichung iiberein- 
stimmen. Die weiteren Glieder in (4. 16’) entstehen als Resultat der Annaherung 
und kénnen vernachlassigt werden in dem Sinne, dass man die Summations- 
grenze m geniigend gross wahlt, um die erwiinschte Genauigkeit zu erhalten. 


Von (4.15) findet man, dass in jedem Falle, das heisst, was immer auch der 
Wert von m ist, 


; (4. 19) 
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_und daher erhalt man fiir die angendherte Lésung von (4. 11) 


m at 
Dares ‘ | 


3 =] Aji a = Ai-1 ; 
Bes ian |S 7-91 a Aa 


j=l it 
DS en 


t=1 


4 

Bbieser Ausdruck bestimmt fiir jedes 2, das die charakteristische Gleichung er- 
 fiillt, die entsprechende angendherte Schwingungs-Kriimmungs-Linie. 

; Lasst man jetzt m —> oo, was augenscheinlich erlaubt ist, da alle Reihen kon- 
_ vergieren, so findet man leicht 

a 

4 

; Zi(y) = A Aare sinh (A1/4 ») — sin (Al/4 g) 

; (4.21) 
y 
4 
2 
y 
: 
; 
: 


1 os + cos Até4 


= | [cosh (42/4 y) — cos (41/4 g)] iP 


sinh 41/4 — sinAl/4 


_ das heisst die seit langem bekannten exakten Lésungen. 


5. Schwingungen des eingespannten Keiles 


Dieses Problem ist im Jahre 1879 durch G. R. KircuHHorr!) exakt gelést 
- worden, indem er von der Differentialgleichung ausging und zeigte, dass die 
Schwingungslinien mit Hilfe von Bessel-Funktionen ausgedriickt werden kénnen. 
In diesem Falle sind sowohl die Massen- wie auch die Steifheitsverteilung 

_ Funktionen der Langskoordinate des Balkens. 


Fig. 4 

Setzt man fiir die Breite 1 
: 1) = 4-3-9, (5. 1) 
my) == 72 a 9), (5. 2) 


1) G. R. Kircunorr, Uber die Transversalschwingungen eines Stabes von verdnderlichem Quer- 
schnitt, Berliner Mber., 1879. 
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so erhalt man von (2. 6) 


3 Pde. — 
Dy: v—»?[ ees a]. (vy > 9) 


Ht) ==78 ne (5. 3) 
d- [5% +@-»]. o<o} 
Nach Einfiihrung der Koordinaten (4. 6), (4. 7) gibt dann (3. 3) 
to oy, ( 
oa p> ¥) 
= a (ole 1 i (° 2 (5. 4) 
K(9, y) = Aion | gel? pa! ; 9 ‘ty 
? (v = 2). (p<y) 


Somit hat nun die Integralgleichung (3. 2’) die Form 


-3| [gale Einav {Soe (r— Z} 29 ay]. (5. 5) 
- LO 


2 


wobei 

3014 ow? 
Eb 
Der Integraloperator auf der rechten Seite transformiert die Potenz y” in den 
Ausdruck: 


A= (5. 6) 


gun gel? gp! 


(F(t F)+ S45) 2-3) ceed) 


Wahit man 


(5. 7) 


A y+ B ys! 


als Anfangsfunktion, so lasst sich das Resultat von m Iterationsschritten in der 


folgenden Form darstellen 
m 


Z(q) = Di {Az p48 + By ptitiay (5. 8) 
q7=1 
Man postuliert jetzt den Fehler zu 

Qa sds 

7h ele 

(2m —1) 2m (2m + 2) (2m + 3) 
pemtsie (a 
+AB 


™ 2m (2m +1) (2m + 35) (2m-+ 4) 


und erhalt nach Anwendung von (5.7) aus dem Koeffizientenvergleich der Po- 
tenzen von 


4;= aT A - ! 
j (27-4 1)1 (27 = 2)) ** B; (@742)(@7—1t (5. 10) 
A,+ . j 4 j Dao 
a 24 6X ati | 
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Die Ausdriicke fiir A; und B; (5.10) in (5.11) eingesetzt, ergeben ein homogenes 
Gleichungssystem in A ‘und B, “~wobei die ausmultiplizierte Determinante fiir die 
picnitriviale Lésung die folgende Form annimmt: 


We 


m j 
—1 
ore (27+1 1(27—2)!2]4 aes aes yl2q! 


2) ll aK 
: ecient (5. 12) 


2) arsait ai eater eee 


Nach steigenden Potenzen geordnet, sollen nun die ersten » Glieder mit der 
Losung von KrrcHHorFF (5. 13) 


1)9 a3 


oe) es a 
& aes ° aie) 


iibereinstimmen. Um dies zu beweisen, ordnet man die Doppelsummen nach 
Potenzen von / und zeigt die Identitat der (m + 1)-ten Potenz von / in den beiden 
Ausdriicken (5. 12) und (5. 13). 


Gj 
E 
4 
A 
7 
: 
. 
a 
E. 
4 
; 
; 
— 
; 
$ 
yf 


m qn 
+1 a. ant? 2)12n! (2-4 2) (2" = 1) — 2] 


. 


; m—1 1 
ie OS area Zl[2(m—K+1)—1]! (2K +1)! (2K —2)!2K 
m—1 n 1 
: se 2 a es (2(w—K 41) +1)! (2 —K 41) 1! @KFA1(2K—1)! (2K +1) 
z ar m qn (—1)" ‘ 
. =D) Gabi mr yi@s yy!’ 
1 (5. 14) 

(2m+4)(2n+1)—2 

(2 + 4)! (2+ 2)! 

yee Pe 
+2) gees ke ERS fo k= 1) on 
: De 
—(2m—-K+i1)+i!(2m—K+1)—1!(2K+2)1(2K+ 1)! 


= ett i . 
— (2 + 3)! (~+-2)! (n+ 2)! 


ZAMP V/5 
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Nach einigen Umformungen nehmen diese Ausdriicke die folgende Form an: 


2n+2 


lara er aleat ss eee 


Setzt man noch m = n+ 1 und beachtet die Beziehungen 


PAs (" "Ke ‘) eek Fame (—1)" (1 ip 


Bey Ne oe cone. 


K=0 


so ergibt sich sofort, dass die geforderte Identitat erfiillt ist. Der Fehler (5. 9) 
verschwindet beim Ubergang m-—> oo, so dass die Lésung ebenfalls die Form 
(5. 13) annimmt. 

Die Kriimmung ergibt sich aus den Beziehungen (3. 3), (5. 8) und (5. 10) zu 


m ' p2i-2 ! m2d-1 ‘ 
ao = 2 Tar ei ee 
mit 
K= ; psa (5. 19) 
| z Dd E 
Durch Wahl der neuen Variablen 
Ger quan 
und der neuen Konstanten 
KA, 31/=As BB, KB) 41142 =]A—B (5. 20) 
erhalt Z(¢) die Form 
2m+1 j 2m+1 a4 1)3 
eee weart® a G Sire 


die mit der Kirchhoffschen Lésung iibereinstimmt. 


Summary * 


The vibration problem of certain dynamic systems with polynomial mass and 
stiffness distributions can be expressed as a Fredholm integral equation with a 
degenerated, symmetric kernel. 

If a starting function is chosen appropriately, the eigenfunction can be ex- 
pressed as a power series. Simple recurrence relations between the coefficients 
of this power series yield the characteristic equation for the eigenvalues with a 
finite number of disposable coefficients. 


This method is applied to a beam and a wedge and leads to the exact solutions. 


(Eingegangen: 21. Februar 1953.) 
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Sur la thermodynamique des processus irréversibles 


(Troisitme mémoire)?) 
Par KyRiLLE Poporr, Sofia?) 
_ Considérons un systéme, adiabatiquement isolé, dont l’état est compleéte- 


ment défini par les variables é,, &, ..., &,. Soient &9, 8, ..., 2 les valeurs de 
ces variables, correspondant a l’état stable du systéme, et posons x; = &; — £9, 


Le a 


¥ ? : ? ° yd ra LJ x 

_ L’expression de l’entropie autour de l’état stable, réduite A ses termes du second 
3 degré, est donnée par 

4 1 , 

AS =~ DY) bia He Xn s avec Sin = Beis @, R= 1,2,..5, 2) 

4 1,k 

L, 

eo > Six %; X%, est une forme quadratique, positivement définie. 

i,k 

_ La théorie des processus irréversibles de la thermodynamique, esquissée par 
_ ONSAGER, repose sur la symétrie du tenseur L;,, figurant dans les relations 
_ phénoménologiques 


x= DT Lin Xe, (1) 
D 


ad a baal 


ou les x; sont les «flux», c’est-a-dire les dérivées des x; par rapport au temps /, 
_ et les «forces» X;,, les dérivées de (—AS) par rapport aux x, respectivement. 
_ Dans ce qui suit, nous montrons: 

| 1° Si les x,;, et par conséquent les g;,, sont d’aprés l’hypothése les seuls 
éléments déterminant l’état du systéme, les L;, doivent étre complétement 
_ déterminés par les g;,. L’intégrale générale du systéme d’équations différen- 
4 tielles (1), qui est d’ordre x, contenant 7 constantes d’intégration ; ces constantes 
4 sont complétement définies par les valeurs initiales des » variables x,;, mais 
_ pour que ces intégrales déterminent des processus irréversibles, il faut qu’on ait 
_ X;(+ co) = 0, ce qui impose nouvelles conditions aux L,; ,. 

. 2° Nous montrons que toutes ces conditions sont remplies par des intégrales 
_ du systéme plus général d’ordre 2 


» » tea wil 


4 


; 0(AS) 
ou Xx; = Sata? (2) 


z 


a*x, 
ae 


dans le cas le plus général, les L;,, ne dépendant ainsi que des g; ,. 


1) Cf. Kyritte Pororr, Sur la thermodynamique des processus irréversibles, 1&€ mémoire: 
_ZAMP 3, 42-51 (1952); 2e mémoire: ZAMP 3, 440-448 (1952); Sur les relations phénoménologiques 
_d@Onsager, C. r. Acad. Sci. Paris 253, N° 18, 648—649 (1952); Sur l’échange de chaleur bar conduction 
aun systéme a un autre, C. r. Acad. Sci. Paris 256, N° 2, 785—786 (1953). 

2) Institut de Mathématiques a l’Université de Sofia. 


68 KyRILLeE Poprorr > ZAMP 


3° Nous appliquons la théorie générale que nous développons a la théorie de 
la conductibilité de la chaleur et a la théorie des phases. 


Etude des relations phénoménologiques d’Onsager 


Considérons le systéme d’équations différentielles 
ay = Ly, X, + Lyp Xn t+ 01° + Lyn Xn, 


Hy = Ly, Xy t+ Lyn Xn +17 + Len Xn, (1) 


@ieviehel-a. 42 e0iejey se) #1 le fo er e116) 8) On etie\ ie) el iol oe, shies, (ote: iP aS us 


avec 


Xe = Sur %1 + Bue %2 + °°° + Gen Xn - 


C’est un systéme linéaire d’ordre 1 dont la solution générale est donnée par 


x,(t) = C,a, ce’ + Cage +--++C,a,e™, 


X(t) = Cy By emt 4 Cy Bo eT fare f Cape ee Fe 


70, 2.0) 6 af eyie, oe, (| ev e's] «6. 0) OS) @ re) @ sm n> eee 6) Sim ee) Wie @ 6 «1 


oti les 7, et les «,, B;,...,A; sont des fonctions des L;,, et des g;, seulement. Les 
conditions initiales x,(fj) = x9, %o(t)) = x8, ..-, %n(%p) = x8 déterminent com= 
plétement les 1 constantes d’intégration C,, C,, C, en fonction des L;,. Mais 
pour que les équations (1) définissent un processus irréversible, il faut qu’on 
ait encore %,(+ co) = 0, x,(+-co) = 0,..., x,(+.0c0) =0, ce qui exige que tous 
les7,soient négatifs. Or, les 7, étant des fonctions des L;,,, cela impose ” nouvelles 
conditions aux L,,. Ainsi, dans le cas de m = 2, ona 


1 
eg ey (Ly3 811 + 2 Lye B12 + Loe 8o8) 


1 = = 2 
ciency Vic Bir + 2 Lie 81a + Lee Seo)” — 4 (Lys Le — Lio) (811 822 — Sia) - 


On constate d’abord que les L;, entrent comme éléments déterminant le 
processus irréversible, contrairement a l’hypothése admise que les &; sont les 
seules variables déterminant l’état du systéme. Il faut de plus que les 7, soient 


tous réels et négatifs pour qu’on ait x,(+ co) = 0, ce qui impose des nouvelles 
conditions a remplir par les L; ,. 


el a a 


OG ee ee ee Re ee eee eee ee ee, ' 
q . : — ere ey 
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Les équations (2) déterminent complétement un processus irréversible 


Considérons maintenant le systéme d’équations différentielles d’ordre 2 


a? x; 


aha X:, @=1,2,...,n) (2) 


_ qu’on peut écrire sous la forme 


d*x, : 
gee 811 %1 + Bie Xe + Bin Xn» 


Ate 


qe 821 %1 + Ben Xe + °° + San Xn» ; 
(2 bis) 


SP b (ele) 616) 6 (6 (6 © (6 (6.0 6 © ole 0 ev ‘© © © 0 106 0 « 6 6 8 6 


ax, 
Vike. = 8ni %1 1 Bue Xa 1 °°" + ban Xn - | 


Ces équations admettent des intégrales de la forme 
iene be ty = P COb a. ky Se, 
ot: les «, B, ..., A et r sont déterminés par le systéme d’équations algébriques 
i) hh Ste P + n=O, 
Ber %=(Boe — 77) B+ °° + Sen = 0; 


le ie: (610 w wget se) Bille ele) 6 B08) 9) (0) 4) 6 101 o lee, 0:16, @, 0) 4: 0110 


Page fas be nae © *)A=0 


homogénes par rapport aux «, f, ..., A. Par conséquent, on a pour déterminer 7 
Véquation algébrique suivante: 


fae" bi fin 
821 822-1" San at (4) 
gn1 gre Enn — " 


qui est de degré 2. La forme quadratique Py Lin Xj X, etant positivement 


_ définie, les m racines (r?) de cette équation sont toutes positives, ce qui fait 
que les m racines de cette équation sont négatives et les autres » racines sont 
~ positives. Soient 7, 72, ---, %, les racines négatives et-%n41) Yn) +++» Yen lesracines 
- positives. Dans le cas Ae n = 3, par exemple, on aura pour déterminer 7? 
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l’équation 
yo — 7* (215 + bes + Ess} 


811 812 §13 


| — | &e1 Gee Ges | = O0- 
831 832 §33 


£11 813 §22 823 811 812 | 
+ + 


st fo 


§31 833 §21 S22, 


£32 §33 


(— AS) étant une forme quadratique positivement définie, les g;; et tous les 
déterminants qui figurent dans (5) sont positifs. On a par conséquent trois 
racines 7? positives, ou bien trois racines 7, 72, 73 négatives et trois autres 
4, 15) %¢ positives. 

Dans le cas général, aux 2 ” racines de (4) correspondent 2 m systémes de 
solutions (x, B,..., x) des équations (3), ot l’on peut supposer toujours les 
Cele 

Ainsi la solution générale du systéme d’équations différentielles (2) est 
donnée par les intégrales de la forme 


%4(t) = Cy me + Cy ge’! +--+ + Ca tne 
ay Cozi tinge doctor Genta ee 
Pour qu’on ait x,(+ co) = 0 (¢ = 1, 2, ..., 2), il faut mettre 


Cia Cris vats Con =0 ’ 


ce qui donne, en posant «, = 1, 
tty (t) = Cy eC, ee 4 = CL em 
a(t) = Cy Bye” + Cy Bye +--+» + Cy Bae™, 3 
6 


Stee es Oa 6k 8 0 ES Ss OS we SC We He Ri ee SS SS ee eS Bee See 


Ces intégrales déterminent complétement les n constantes d’intégration 
Cy, Cy, ..., C, par les valeurs initiales x,(f.) = x°. 

Avec ces valeurs des intégrales x;(¢), on obtient, aprés des réductions faciles, 
pour les X; des expressions de la forme 


X= Cie A+ Ce A, +-++4+ Ce A,, 
xo Cy en! By of Gs et! B, + pace As G. ent a 


Ps <, s (7) 
X, = Ce" N, + Cie" N+ --- +C, ef N, | 


mes | 


SR ee my ee a See eee a re 
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_ otles A, B,..., N comme les «, B, ..., A ne dépendent que des g;,. On a de plus 


[ 


x(t) = C, ee"! 7, Weds Gwe a te Ne +-:-+C,e" 7, Tins (8) 


L’élimination des C,, e”’, ...,C, e’* des équations (7) et (8) conduit aux 
relations phénoménologiques 


Kee P59 da Me ae Tu Ne 

pep ne, Cams ae 

ie Bie Be B, |=0, 

Xe N, Nz NN, 
qu’on peut écrire sous la forme 

ye = Diy Xy + Leg Xyt--s + Lin Xn, (= 1, 2,...50) (9) 

et ot les L;, comme les 7, «, f, ...,A sont des constantes ne dépendant que 
des g;x- 


Etude du cas général ou l’on garde tous les termes 
dans le développement de AS 


Dans le cas général ou l’on garde tous les termes dans le développement de 
AS, ilest essentiel de montrer que l’on am intégrales x,(¢), contenant 1 constantes, 
permettant de choisir arbitrairement les valeurs initiales x? et s’annulant pour 
t= + co. En suivant la méthode d’intégration de POINCARE, que nous avons 
employée dans notre second mémoire cité, on peut ramener les équations (2), 
ou cette fois on retient tous les termes dans le développement de (4S), a un 
systéme d’équations différentielles d’ordre 2 » de la forme 


du, ee duy Peer du, a dUnyr een 
Uy + + Yq Un + °°: Yet apes Pog Uns, tex (10) 
— dan = dt | 
Yon Un + tee 


ov les points (...) remplacent des termes de second degré et de degré supérieur 
4 deux, ot les u, sont des fonctions linéaires des x; et des x; et ott les 7; sont les 
racines de l’équation algébrique (4). Aux racines négatives 7, 72, 7, correspon- 
dent ici intégrales x,(t), fonctions holomorphes de ent eft ||. em? s’annulant 
pour t = + oo et contenant  constantes d’intégration et qu’on peut développer 
suivant les puissances de e”’, e”’, ..., e’n!, Avec ces valeurs des x;(¢) on forme 
les X, = —6(4S)/5X,, qui deviennent ainsi des fonctions holomorphes des att, 
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L’élimination des e”* de x{(t) et des X,, XQ, ..., X, conduit ainsi aux relations 
phénoménologiques. Dans ce cas pourtant les L;;, sont des fonctions de ¢ aussi, 
dont les valeurs principales coincident avec les valeurs du paragraphe précédent. 
Pour plus de détails, nous renvoyons au tome III du Traité d’analyse de E. Pic- 
CARD et a notre second mémoire cité. 

Aprés ces généralités, considérons quelques problémes physiques concrets. 


Transférement de la chaleur d’un systeme a un autre 


En appliquant la méthode ici développée nous allons montrer que l’hypo- 
thése de FourIER concernant le transport de la chaleur d’un systéme a un 
autre découle des équations (2). 

Considérons un systéme de deux plaques cylindriques, mises en contact 
thermique par leurs faces et adiabatiquement isolées du reste par une couche 
imperméable a la chaleur. Dans ce qui suit, nous négligeons le travail de la 
dilatation thermique. 

Soit m la masse, par unité de surface de contact, 7, la température absolue, 
c, la chaleur spécifique de la premiére plaque et mp, 73, c, les données correspon- 
dant a la seconde; 7, la température du systéme aprés |’établissement de 
l’équilibre thermique. Nous admettons que l’épaisseur des plaques est assez 
faible pour que la température de chacune reste uniforme pendant le processus. 
En désignant par q, la chaleur 4 emprunter, par l’unité de surface, au moment ¢ 
par la premiére plaque de la seconde pour arriver a |’équilibre thermique et par 
q2 la quantité correspondant 4 la seconde, on aura g, + gq = 0, ot 


Gy = M,C (T, — To) = my Cy AT, Go = Me Cy (1, — Th) = Mz Cy AT, , (11) 
d’ot: 1’on tire 
My Cy Ty + my Cy Tz = Ty (my cy + me Cp) (12) 
et 


My Cr+ My Cy 
My, Cy + Mz Co 


n= (Z,—-Z}). (13) 


On aura ainsi pour AS, correspondant A un moment donné ¢ et A l’unité de 
surface de contact thermique, en ne retenant que des termes du second degré, 


ce qu’on peut écrire, eu égard a la formule (13), 


Aes oO ee 
My Cy > Mr Co Ree 


<3, ev 
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d’ot l’on tire 


0(AS) 2 MC + Mz Cy 
0q, ils My 6) -M> Cy 


%=4q,, avec a>od, 


et, par conséquent, on aura ici l’équation unique 


dq 
ae Ns (14) 


Edont Vintégrale est de la forme g, = Ce". 
En posant cette valeur de qg, dans (14), on aura, a étant positif, 


pst) les a. 


la valeur — Va correspondant seulement au probléme physique, g, devant 
tendre vers zéro pour ¢t = + co. 


4 
- 
7 


; On a ainsi 
m=Ce=Cel", ga=—Yacel*——Vag, (15) 
G 
“avec C = g{—la valeur de g, pour un moment donné ¢ = ép. 
_ Avec la valeur de g, de (13) et de a, on obtient l’expression suivante du 
; «flux» de chaleur, correspondant a l’hypothése de FouRIER 

1 2 My Cy > Mo Co 

A= <5 OS Gen. (16) 


Hc alte in nN Dat be ~~ 


Jusqu’ici nous avons supposé qu’a chaque moment la température des 
plaques est uniforme. Mais quoique trés minces, la température des plaques, 
_ qui est une fonction continue le long des axes des plaques, différera sur les deux 
‘faces de chacune des plaques, et on aura des valeurs différentes de l’entropie 
“suivant qu’on forme l’expression de l’entropie avec l’une ou avec l'autre de ces 
-températures dont la différence est petite de l’ordre de |’épaisseur des plaques. 
_La valeur de l’entropie de chaque plaque correspondra a une température 
“moyenne, dépendant de la conductibilité des plaques. Si l’on désigne par J, et 
eT, les températures extrémes sur les deux plaques, 7) la température a |’état 
Pde stabilité et en posant AT, =7,—7), AT, = T,— Th, les températures 
-moyennes A mettre dans l’expression de 4S seront 


- 


Renn AT Tet A Alas 


< 


ot A, et A, dépendent de la conductibilité des plaques et différent trés peu de 
‘Vunité. Ainsi on aura 


q = m4, 4T,, Jo = Mz lglg ATz, 1 +%=9. (17) 
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Avec ces valeurs de q; et g, on aura 


= My Cy Ay + Mg Cy Ay cp Tt. 1 UAT 
a= My Cy Ay + Mg Cy Ag (1 2)» my alien (18) 
Oo ne eis yh pea My, Cy + My Cy 
me Ms Cy’ Ge a asta 1 my Cy + My Cy 
et en se bornant aux termes de second degré 
AS Nn a Yo 24s eS My, Cy + Mz ly Gi 
t,o Aad Ty + A, AT, Te . IM, Cy + MaDe Pa 
d’oti 
A(—AS) 2 (m+ Me) aq 
aq, = My Cy - My Cy TS i 


En suivant les calculs de la page précédente, on aura définitivement 


aby ae MGA My Cghg TT, —-T, 
a= My Cy + Mz Ce My Cy Ay + My Cy Ay Loy 


Conduction de la chaleur dans une barre homogéne 


Considérons une barre cylindrique de densité 9 et de chaleur spécifique c, 
mise dans un étui imperméable a la chaleur, éventuellement ouvert méme de 
deux cdétés. Nous admettons que la température de la barre est une fonction 
continue de l’abscisse x, comptée sur l’axe de la barre, et admet des dérivées. 
Considérons trois sections de la barre, correspondant 4 x —dx, x et Aa x+dx 
et soit f l’aire de ces sections que nous supposons égale 4 «7. On a ainsi deux 
‘tranches limitrophes de volume / dx et de masse @ f dx. 

Considérons le systéme des deux tranches. Soit g (x — dx) le flux de chaleur 
pénétrant le systéme par la section (x — dx) et g (x + dx) le flux de chaleur 
sortant par la section (x + dx), rapportés a l’unité de surface. 

Désignons par ds l’accroissement de l’entropie de l’unité de masse du 
systéme des deux tranches pendant l’intervalle de temps di. On aura pour 
l’accroissement de l’entropie du systéme, rapporté a l’unité de temps 


ds 
2 of ax re 
Désignons d’autre part par du l’accroissement de I’énergie interne, rapporté 


a l’unité de masse pendant l’intervalle de temps dé. On aura pour l’accroissement 
de l’énergie interne du systéme pendant l’unité de temps 


d : : j 
2ofdx F=f [9 (x — dx) —q (x+ dx] = —2 fax OO. 
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Ainsi, en négligeant le travail de la dilatation thermique, on aura 


ds _ dq(x) 
Te CO er eer rome 
- ou bien 
2a Mag) q_, aT(#) 
Se paaeete ae Fe dx TE dx 
~ Pourtant ona 
; zs 
F T(x + dx) —~ T(x — dx) =2 97) gy 
4 
4 Or, d’aprés la formule (13), cette quantité est proportionnelle a la quantité 
4 de chaleur g dans un processus adiabatique. En posant 
4 aT 
: aa. te 
4 
-  étant le coefficient de proportionnalité figurant dans (13) ou (18), on aura 
ds d(q/T) qq 
; Cat dx T? 
Ainsi l’accroissement de l’entropie se trouve décomposé en deux parties: 
1° a la divergence de l’entropie 
; 
E wg) 
4 dx 
j et 2° 4 une partie 
4 L 
Ao =— 4599 
d correspondant a 
3 1 lw 
pga: =o q? 


vy ees 


_ et provenant d’une source, due a un processus adiabatique irréversible, consi- 
_ déré dans le paragraphe précédent. La dérivée par rapport a q de cette partie 
— étant (—y g/T?), on aura pour déterminer g, d’aprés la théorie développée ici 
; a*q ue 


"de = b q avec b = care 


_et qui conduit, d’aprés le paragraphe précédent, a 


; § sie eg 
g--Voq--V/5 Nida 
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C’est la formule du flux de chaleur, d’aprés l’hypothése de FourrErR. Dans son 
ouvrage Thermodynamics of Irreversible Processes, M. DE GROOT reproduit le 
résultat analogue, di 4 M. PricocIne et obtenu par l’application directe des 
relations phénoménologiques d’OnsAGER, dans leur forme primitive ou figurent 
les coefficients L,;,, dont la signification est, comme nous l’avons montré, a 
établir pour qu’on ait ¢(+ co) = 0. 


Applications 4 la théorie des phases 


Considérons un récipient de parois imperméables a la chaleur et absolument 
rigides, rempli d’un liquide (I) et de ses vapeurs (II) a des températures respec- 
tives T, et Ty. Soient M, V, U la masse totale, le volume et l’énergie totale de 
la substance, remplissant le récipient. En désignant par M; (2 = 1, 2) les masses 
respectives des deux composants, par v,, “, et s; le volume, 1l’énergie et l’en- 
tropie respectifs, rapportés a l’unité de masse des composants, on aura 


M,+M,=M, M,y,4+M,x2,=V;, M,44+Miuw=U, | 
(20) 
Sa irs ios | 


Choisissons l’énergie et le volume comme variables indépendantes déter- 
minant l’état des composants. Les relations (20) déterminent Mg, v,, wu, en fonc- 
tion de M,, v,, uv, D’autre part s;, T;, p;, étant des fonctions de w,, v;, devien- 
nent des fonctions de Mj, v,, u,. En définitive, l’entropie totale S du systéme sera 
une fonction déterminée de M,, v, et u,. En désignant par 9, v9, M9 les valeurs 
de ces variables a l’état d’équilibre thermodynamique, on aura, en posant 


y-v=v, 4—vW2=4, M,—-M,=m™, 


AS = Sty % My) ~ Slut, of, MY = (52°), + (5>),0+ (sar), ™ 
to) 


Ou, Ov, 


2 


1 7/a2s 02S 02S 25 
ee cee = 2 eee is 
May ( Our i: I Ge alk ee leer ), ee ( Ou? |) BS 


02S 02S \ 
ah, ) a Se 2 Toe 
on aie) 84 Gade: 


Les dérivées partielles étant prises pour les valeurs u9, v9, M ? des variables 
indépendantes, correspondant au maximum de S, on aura 


(35), (35), = (air. -° a 


et en ne retenant que les termes du second degré en u, v, m 


1 DENS 02S 2 
A Ie ae a) 249 [ a 02S 02S 
2 Ouz ue + Ou, 00, Uv + 2 peak um + ee ) y2 


2S Oh Sy mage 
nie one OM, a aa Gree m?| 
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Pour le calcul des dérivées secondes, nous aurons recours a l’expression de 
la différentielle seconde 


DENS 07S GES) 


Aor a 
as = ai Oe a en 


Les équations différentielles, correspondant aux équations (2), sont dans 
le cas 


du as A 2S 02S 
aie ( dur Jo Ou, Ov, y vo ( du, OM, ), ae 

dv 02S 02S 02S 

dt? ( Ou, 00; ), #1 ( Ov? iL X! ( dv, OM, bi ie (22) 
dm 02S 2S 2S 

dt? ( du, OM, ), a ( dv, OM, ih : ( FIVE ), he 


dont la solution qui nous intéresse dépend des racines négatives de l’équation 
en 7 suivante: 


02S oa ye 02S 
e a i; a Sala 00, ), du, 0M, Jo 
. 02S 02S i 02S = 
- ( Ou, OV, ), ( Ovt i o ( dv, 0M, 1) =0. (23) 
02S 02S - ( 02S es 
aj en ah ~ \du, 0M, Jo 0M? }o 


; 
, 
7 
3 
, 
4 
: 
; 


—(AS) étant une forme positivement définie, trois des racines de cette équation 
en r sont négatives et les trois autres positives. Les constantes d’intégration 
dans l’intégrale générale du systéme (22), correspondant aux racines positives, 
doivent étre posées nulles. 

Les résultats qu’on obtient ainsi correspondent au probléme physique tant 
que la masse d’un des composants, donnée par les équations (21), n’est pas 
négative ou plus grande que M. 

Pour déterminer les valeurs des variables indépendantes correspondant a 
V’état d’équilibre thermodynamique nous suivons la méthode de Gisss, telle 
qu’elle a été exposée par Max PLANCK). 

Des équations (20), on obtient 


Vere 


dM, = —aM,, M, dv, = —M, dv, — (v, — U,) 4M, , 


SS Rs 8 bd ee Bl 


M, ditty = —M, du, — (uy — %) 4M, 


1) Max PLANCK, Thermodynamik, 9¢ éd. (W. de Gruyter, Berlin et Leipzig 1930). 
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et : 
dS = M, ds, + M, ds, + s,4M,+s,dM,, (dM = dm) 

= M, ds, + My, ds, + (s, — 5.) dM 


SM: du, ap dv, M, du, + py dg + (s,—s) dm ost 
eae (= es z) du, + M, ( a Bt) dv, 
+ (s: Ss sae 7 kee 2) ) dm . 

Ainsi les équations (21) définissant l'état d’équilibre, deviennent 
Bm (h-H-o, Bom@-B)-0. | 
0S _ gt + (4 —%) 9 Ge) 
Ce ag aie iz 

d’ot l’on tire les conditions de GIBBS 
T= len Pe Pa Pe: (27) 


D’autre part, u et v ayant été choisis comme variables indépendantes, on 
aura 


Ty = filtty, dy), Ty = f(a, V2) = (tty, Vy, My) , 
Pr = Pil, 1), Pa = P2(Ue, Vo) = PM, %, M4) , 
Sy = Yi, Vy)» Sp = Po(My, Vg) = (tty, Vy, My) 


et les trois équations déterminent les trois variables v,, u,, M, qui correspondent 
a V’état d’équilibre thermodynamique. I] est 4 remarquer que les équations (26) 
contiennent s et w par la différence de leur valeurs a l'état (I) et al’état (II). 

En suivant Max PLANCK, on peut donner une autre forme a la derniére des 
équations (26). La différence (s, — s,) a l'état d’équilibre thermodynamique ne 
dépendant que de I’état (I) et de l'état (II) est indépendante du chemin suivi 
dans le passage de l'état (I) a l'état (II). On peut trouver cette différence en 
suivant le passage sur l’isotherme 7. On aura ainsi, s, et s, se rapportant al’unité 
de masse de la méme substance et 4 la température Jy, 


- d iar 
uU + LU U, — Us 
Sy Sg = [as =/ — ae => mare — ao T i dv A 
0 
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Avec cela, en posant T, = 7), la derniére des équations (26) devient 


[ba a 0, [pdm pt 1- %)- 


Il nous reste former les dérivées secondes, figurant dans les équations (22), 
au moyen des quantités mesurables directement. 


B De at on obtient, eu égard aux relations (21) 
@S = + dT;) du, + a re et aT do, eo 7 (dp — dp) dv, 
we _ du — du,+ he (dv, — dvs) (vy — 2») is 
+ A ds, T, a T, 
r Uy — Up (vy 2) aT,| dM, 


= 3 (—aT, + aT) du, + “ye (j aT, + 4T,) doy + (dp, — dp,) dv, 


i es = sve 
if e Us i (vy; — Ve) oi. = cae ae dp,| dM, . 


j 
4 
; 
r 
| 


Les 4, 01, U2, Vv, et M, qui figurent dans cette expression de d?S correspondent 
r a l’état d’équilibre. Les différentielles d7;, dT, du, dv,, dp,, dp, d’autre part 
4 doivent étre exprimées en fonction de du, dv, et de dM, = dm. 

Nous avons eu d’abord 


E 
a 
4 M, dv, = —M, dv, = (v, — v,)dm, avec M,=M—M,, 


wet v ayant été choisis comme variables indépendantes, on aura T = T(u, v). 
| Cette relation détermine ~ comme fonction de T et de v, soit u = u(v, T), 
‘satisfaisant identiquement l’équation T= 7[u(v, T), v], d’ou I’on tire ae 


en prenant les dérivées des deux membres suivant T 
3 ‘ 
E- 1 Ol, v} du(v,T)° - OT(u,v) 1 1 

: = 


ee oi Ob duu, ToT e~ 


v 


En différentiant la méme identité suivant v, on obtient 


OT(u,v) _ _, * OT (u,v) | Outv,T) __ i du(v, T) 


a du dv fo 08 


_ D’autre part, on a d’aprés le premier principe de la thermodynamique 


dg = du(v, T) + plv, T) dv = OD) ar 4 [MOD + piv, T)] dv 


T) 
=c,dIT +adv, avec a= el t) + by, 10 
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et par conséquent 


OT (u,v 1 
“ L = — re (a—?). 
Définitivement, on aura 
OT (u, v) OT (u, v) a p ce 
I = a du + — d Z dv + u 
et par conséquent 
Pi eee ELOY Pe ae yeh 
vl v 
— 1 
Fd Gia pe Se LN SP nD 
Cy2 Cy2 
a—p, M, <p 1 M, 
Cy : M, an Cue My, oe 


Pour obtenir dp en fonction de du et dv, considérons # comme fonction de 
u et de v. On aura ainsi 


ap(u, v) = aS ») du + = ”) dv : 


D’autre part, on a d’aprés le premier principe de la thermodynamique 
dq = du(p, T) + pdv(p, T) =kdp+c,dT 
ou l’on a 


vP.T) dup, T) , , Ov(6, 7) 
Op noe tp 2 yee rv 


En portant dans l’expression de #(w, v) en fonction de wu, v les valeurs de 
u et v en fonction de p et T, on aura 


p = p(u,v) = plu(p, T), v(p, T)], 


d’ot l’on obtient 


ie Op(u, v)  du(p, T) oh Op(u, v). dv(p, T) 
Ou Op ou Fp pees 


ie Op(u, v) ; du(p, T) Op(u,v)  dv(p, T) 
Ou oT oD 7 Son 


Vol. V, 1954 Sur la thermodynamique des processus irréversibles 81 


_ D’ici, en tenant compte des expressions de & et de c, et en posant 


ie dv(p, T) Ou(p, T) 
3 A=k—JF Pope 
- on obtient 
a Op(u, v) Se du(p, T) —, 0ulp,T) Op(u, v Ov(p, T 
eigee tt UST PS pr ie = 
dott 
ap = ae v) PRS ae v) heed eee, vee Lee dv. 


On aura ainsi 


Beg ME UIE agg (se BLOM TOT gy 


a bts 


et en tenant compte des expressions de du, et dv, 


M, | 0v(p, T)/oT = , 
Fe EEE) ga 5 (Se Pe DPT 


1 Ou(p, T 
— aa; {la — + Pol Vp) | ies ey (vy — %»)| | dm, 


les indices 1 et 2 se rapportant aux phases (I) et (II) a 1’état d’équilibre. 

De l’équation dq = du + # dv, on obtient pour la quantité de chaleur (7, — q,) 
nécessaire pour faire passer l’unité de masse de la substance de température T, 
et sous pression f, de l’état I a ]’état II, 


aE ey Me NTE nee a ee ne ey eT eh ere eee eee 
er : 5, 7. 


Ia — a = Uy — Uy + Po (V1 — %) - 


4 Ces valeurs des dT et des dp introduites dans l’expression (28) de d?S con- 
4 duisent, aprés réductions faciles, 4 une expression de la forme 


a2S = A du? +2 Bdu,dv,+2C du,dm+ Ddvi+2E dvudm+F dm? 


_ dont les coefficients sont les dérivées secondes respectives figurant dans les se- 
- conds membres du systéme d’équations différentielles (22), Nous nous dispensons 
_ a reproduire ici les expressions de ces coefficients. Nous nous contentons de 
-remarquer qu’ils sont formés tous avec des quantités qu’on peut obtenir par 
_ des mesures directes. On aura ainsi le systéme d’équations linéaires d’ordre 6 


da?u dv 
ae a Bu Cm, wii eee Di Em, 
d?m 
ee OE Pe, 


_ ZAMP V/6 
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dont l’intégration ne présente pas de difficultés. Les intégrales cherchées corres- 


pondent aux trois racines négatives de l’équation algébrique de degré 6 en 7. 


Zusammenfassung 


In zwei friiheren Abhandlungen haben wir uns mit den phanomenologischen 


Relationen * 
#, = S Lin Ae CR ene nes) (1) 
Rk 


beschaftigt, wo x4 die « Stréme» und X, die « Krafte», das heisst die Ableitungen 
von ¥; in bezug auf die Zeit ¢ bzw. die partiellen Ableitungen von (— AS) in bezug 
auf x, bedeuten und wo (4S) der Ausdruck der Entropie in der Nahe des Stabili- 
tatszustandes eines adiabatischen Systems ist. 

Fiir ~ = 2 haben wir gezeigt, dass die phdnomenologische Relation (1) ein 


System erster Integrale von dem Gleichungssystem 
d*x; : 
ae j (iia 1, 27a) (2) 
ist. 
Hier betrachten wir den Fall, wo » irgendeine ganze positive Zahl bedeutet, 
und zeigen in dem Falle, wo (4S) der reduzierte Ausdruck 


AS = =— Sie tite) beet 
t,k 


der Entropie bedeutet, dass die Lésung des Systems (1), das ein lineares System 
n-ter Ordnung ist, nur » Integrationskonstanten enthalt, welche die Anfangs- 
werte von #,(¢) willktirlich zu wahlen gestatten. Nun miissen aber die Werte von 


n 
#,(t) = Sheers, 
j=1 


fiir ¢= +00, gleich Null sein, und folglich miissen alle x Werte von 7; negativ 
sein. Da die 7; Funktionen der L;, sind, so erscheinen die L;, als Elemente, die 
den Charakter des irreversiblen Prozesses bestimmen und ” neue Bedingungen 
ergeben, was der Grundannahme widerspricht, dass die x; die einzigen unabhan- 
gigen Variablen sind, die den Zustand des Systems bestimmen. 

Im Gegensatz jedoch zu dem ist die Lésung des Gleichungssystems (2 ~)-ter 
Ordnung (2) von der Gestalt: 


2n 
x, (t) = DiC; UE erst 5 (2 = lie a 3; ais ay n) (3) 
jai 


wo nun die v7 die 2” Wurzeln der sékularen Gleichung (2 )-ten Grades 


2 
Same §12 CO Siler | 
2 
§21 Sigh 2 S2n 2G 
ihe Pasty siausentetos shodehonetece Renner: Riera 
le g wer — 2 
onl n2 Snn 


sind und nur von den g;,, abhangen. Da ¥ g;, *; ¥, eine positiv definite quadrati- 
sche Form ist, sind alle Wurzeln reell: ” positiv und v negativ. Setzt man nun alle 
C;, die den positiven Wurzeln entsprechen, gleich Null, so bekommt man eine 
Loésung *,(t), die nur  Integrationskonstanten enthalt, welche gestatten, die 
Anfangswerte #° willkiirlich zu wahlen, ausserdem ist hier x, (+4) = 0. 


PA ee 
7 


Pee VAD eae ee 


ae ey) Ne 


AVA pe 
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Die Elimination der C;, e’s¢ aus x,(¢) und aus den mit (3) gebildeten Werten von 
PX; (k = 1, 2, ..., n) fiihrt zu den phanomenologischen Relationen (1), wobei die 
L;;, konstante Gréssen sind, die nur von den &;~ abhangen. 

Es wird weiter gezeigt, dass auch in dem Falle, wo (AS) alle Potenzen der »*, 
enthalt, die phanomenologischen Relationen (1) eine Folge von (2) sind, wo nun 
die L;,, auch von ¢ abhangige Glieder enthalten. 

Von diesem Standpunkt aus sind die Fourier-Formeln fiir den Warmeaus- 
tausch abgeleitet und die Gleichungen (2) auf die Phasentheorie (zwei Kompo- 
nenten) angewendet. 


Regu le 20 février 1953.) 
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_ Note on Two-Dimensional Indentation Problems of a Non-Isotropic 


Semi-Infinite Elastic Medium 


By BIBHUTIBHUSAN SEN, Chandernagore, India 


Introduction 


Rectangular indentation on the boundary of a semi-infinite isotropic solid 


_ was considered by Sapowsky (1928), and indentations of various shapes for such 
_ solids were discussed by SNEDDON (1951) who obtained the result in each case 
_ by solving dual integral equations. It was, however, shown by the author (SEN, 


1946) that in cases of rigid punches on the straight edge of an isotropic semi- 


_ infinite plate, the solutions could be easily obtained by proper choice of complex 


functions in terms of which the stresses and displacements were expressed. In 
this paper, this particular method has been applied to solve the problem of an 
approximately circular indentation on the straight boundary of a semi-infinite 
plate having a simple form of aeolotropy. The extension of the method to other 


_ forms of indentations and more general types of aeolotropy is not difficult. 


1. Stress-Strain Relations and the Expression for the Stress Functioh 


Let the origin be taken on the straight edge along which the axis of # is taken. 
The axis of y is drawn into the plate perpendicular to the edge. In the material 


_ considered we take 


E,, E, = Youne’s moduli in w and y directions respectively, 


o, = Ratio of the contraction parallel to #-axis to the extension 
parallel to y-axis, 


thea 
0, = Ratio of the contraction parallel to y-axis to the extension Ved) 
. parallel to x-axis, 
“ = Modulus of rigidity associated with the directions of ¥ and y. 
We further assume that 
te ey” 1.2 
EB, IDA : (A. 2) 


say, a relation which holds good in a material like wood. 
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The stress-strain relations are 

ban = EFF —RYD, 
wy = REX + Ex* VD, (1. 3) 
Cny = eK. 


Equations of equilibrium in the absence of body forces are satisfied if we write { 


é 


— 074 ron 074 nt ee ee 07x 1 4 
4M yh ee igh Cee aye i 
From the relations (1. 3), (1.4), and the compatibility equation 
OFe é O7G,,  * O7lry (1. 5) 
Ox? Oy? Ox OY 
' we obtain the equation satisfied by x as ; 
Oty Oty O4y 
2) ! 5 6 
Syke’ ro rt One ade Be 
in which 
2M = (E,E,)'?(u-t- 2k), n=ey, e=(E, Ez¥*. (1.7) 
Hence we can write 
4 = Re[K(*¥ + ta) + Hl +2 By], (1. 8) 
where 
Re = real part, i=/—1, | (1. 9) 
a= (M+ (M?—1)1/2]¥2) B= [M — (M2 — 1) 42] 1/2, | 
Finally putting 
Y= as and y= Pe, (1. 10) 
we get 
~= Re[K(¥ +1719) + Rl¥ + ty, *)]. (1. 11) 


By using the relations (1. 4) expressions for the stresses can now be deduced 
from the above result?). 
2. Method of Solution 


In the particular type of problems under consideration, we shall have #y = 0 
when y = 0. Hence assuming 


~= Rely, F(a +iyy) -—wyF*+iny)), (2. 1) 
we can write (cf. SEN, 1939) 

Re[t 71 v2 {F" (4) — F(a) }), | 

= Re[y, F’(4) — 1 F’(@)], (2. 2) 
4% = Rely y8 F"(e:) — v2 YE F"(a)] , 


® 
I 


= 
“S 
| 


1) For materials like wood which we have considered here, y and ¥, will be real positive 
quantities. 


i aia ae 


ee a ae ae Vo oe aN 


+ ke 


| ehh. | 


aw 


- 
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where dashes denote differentiation with respect to arguments, and 
&=*x+iny and z~z=—=*+ip-y. (2353) 


Neglecting the rigid body displacements, we have the components of displace- 
ments in the w and y directions as 


u = Re[—(yj Ey* + R) yp FP’ (a) + (v3 Ep? + k) y, F’(@)] 
v = Re[t (yy Yo) {(R v8 + Ea?) yi F’ (eq) — (Rk v3 + Eg?) v3 F(z) }). 


For indentation problems, boundary conditions to be satisfied along y = 0 are 


(2. 4) 


#y=0, Py=0 (l*|>a), v= v(x) (|¥|<a). (2. 5) 
Moreover stresses should vanish at an infinite distance from the boundary. Put- 


ting z= ¥+7y, it is found that the problem reduces to the determination of an 
analytic function F’(z) of z such that 


HH (Rei F'(e) yo = %0(4) |x|<a (2. 6) 
V1 v2 E 
and 
(Ye — ¥1) [ReF”(z)]y-0 = 9, |4| >a. (2. 7) 
If the forms of F’(z) and F”(z) are known, corresponding functions of z, and z, 
can be written down by mere substitution, and the expressions for the compo- 
nents of stress and displacement can be easily deduced. 


3. Approximately Circular Indentation by a Rigid Punch on the Boundary 
et 


(Y2— 41) F(z) = Sa ee [Va? — 22 +72], (3. 1) 
Then 
(VP)y=o = (Ya — 1) [ReF"(z)]y=0 | 
= oe fal <a | (3. 2) 
=0 |v|>a. 


The total load on the pressed portion is of magnitude 


f [Vem R ae? : (3. 3} 


a 


By integration, we have from (3.1). 


‘7 a AN pe ae 
7) ia aay [2 Va? = 2 + asin 1S tial + const , 
so that 
vi 
Uo(#) = ———= (Re{1 F’(2) }y- 
Cogent pa (3. 4) 
Bose Pyre Va 


Sy eS #2+C, when|¥|<a. 
1 V2 2 
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Assuming that v)(*) = 0 when * = + a we shall have 


P : Yat V2 
uw a® V1 Yo Ee 


v(*) = (a* — x4), (3. 5) 


For the value of F”(z) assumed in (3.1) all stress components will tend to zero as 
y > co. The boundary curve of the indentation will be circular since the curvature 


2 
- =— cor =a constant quantity. 
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Zusammenfassung 


Es wird ein Lésungsverfahren besprochen fiir das Problem der elastischen, 
langs ihres Randes durch einen starren Stempel belasteten Halbebene. Die 
Methode wird fiir einen einfachen Fall von Anisotropie am Beispiel eines Stem- 
pels mit annahernd kreisbogenférmigem Rand illustriert. 


(Received: August 17, 1953.) 


Ein Oszillationssatz fiir algebraische Eigenwertprobleme 


Von FRANK WILLIAM SINDEN, Ziirich?) 


Die vorliegende Mitteilung befasst sich mit dem Problem, die Oszillationssatze, 
wie sie fiir Differentialgleichungs-Eigenwertaufgaben gelten, sinngemass auf ge- 
wisse algebraische Naherungsprobleme zu tibertragen. Die Naherung wird nach 
Art der Differenzenrechnung vorgenommen. Die Komponenten des Naherungs- 
vektors sollen also eine Tabelle der entsprechenden Funktion darstellen. An 
Stelle von Aussagen tiber die Anzahl der Nullstellen einer Eigenfunktion treten 
Aussagen tiber die Anzahl der Zeichenwechsel in der Folge der Komponenten 
eines Eigenvektors. 

Definition. Eine Matrix A heisst variationsvermindernd?), wenn fiir jeden 
Vektor x der transformierte Vektor » = A x immer weniger oder gleich viele 
Zeichenwechsel wie x selbst aufweist. 

Hauptsatz. Man betrachte das Eigenwertproblem 


le DIGS (1) 


wo D eine positive Diagonalmatrix ist (alle d;;> 0) und wo A symmetrisch und 


1) Institut fiir angewandte Mathematik der ETH. 


2) I. ScHoeNBERG, Uber variationsvermindernde lineare Transformationen, Math. Z. 32, 321 
(1930). 


| 


q 
3 
4 
; 
4 
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positiv definit ist. A sei noch variationsvermindernd, und alle ihre kodiagonalen 
Elemente @;, ;,1, @;,;-1 seien verschieden von Null. Es gilt dann fiir die charakte- 


ristischen Zahlen: 
Cee lg ee ih 0 


Sind die Eigenvektoren nach der entsprechenden Reihenfolge geordnet, so hat 
der k-te Eigenvektor genau k — 1 Zeichenwechsel. 

Zusatz. Ist A variationsvermehrend und gelten alle tibrigen Voraussetzungen 
des Satzes, so hat der k-te Eigenvektor genau ” — k Zeichenwechsel. 

Da die Eigenschaft «variationsvermindernd» im Grunde genommen allge- 
meiner Art ist, scheint die Vermutung plausibel, dass dieser Satz auch auf Inte- 
gralgleichungen verallgemeinert werden kénnte. Eine solche Verallgemeinerung 
wird jedoch in dieser Arbeit nicht unternommen. 

Mit Hilfe der Tatsache, dass die Folge 


eis Gx (Vs — V3)y Fe (Va — Va), --+) On (Yn — Vin = Cais Ian 


(samtliche a; > 0) immer mehr Zeichenwechsel als die Folge 91, Vo, .--, Vn auf- 
weist!), kann man durch Iteration und durch Wahl von geeigneten vorzuschrei- 
benden «Randbedingungen» allgemeinere Differenzenoperatoren konstruieren, 
welche die Voraussetzungen des Zusatzes erfiillen. Als Beispiel betrachte man das 
Eigenwertproblem des transversal schwingenden, inhomogenen Stabes. Im Inne- 
ren eines Grundintervalls (a, b) gelte 


py)’ =Aqy,. Pl4)>0, g(#)>0, 


und an jedem Ende seien genau zwei von den vier Funktionen y, vy’, py”, 
(p y”)’ gleich Null. Ist dieses Problem auch noch selbstadjungiert und positiv 
definit, so kann man es nach der eben angedeuteten Methode durch ein Diffe- 
renzeneigenwertproblem approximieren, welches ein oszillatorisches Eigenvek- 
torensystem besitzt. 

Ist eine Matrix A variationsvermindernd, so ist die Matrix A*, die man aus A 
durch Multiplikation aller geradzahligen Zeilen und aller geradzahligen Spalten 
mit (—1) erhalt, variationsvermehrend. A und A* besitzen dieselben charakteri- 
stischen Zahlen. Diese « Schachbrett-Transformation » kann zum Beweis des Zu- 
satzes und zur Herstellung weiterer Beispiele benutzt werden. 

Zur numerischen Berechnung von Eigenwerten stehen verschiedene « Ein- 
schliessungssatze» zur Verfiigung. Diese geben Schranken an, die mindestens 
einen Eigenwert eingrenzen. Es wird jedoch nichts dariiber ausgesagt, der wie- 
vielte Eigenwert durch diese Schranken eingegrenzt ist. Der folgende Satz, der 
eine Modifikation des Collatzschen Einschliessungssatzes?) darstellt, liefert im 
Gegensatz dazu Schranken fiir einen bestimmten, durch seine Nummer gekenn- 
zeichneten Eigenwert. 

Fiir das Eigenwertproblem (1) seien die Voraussetzungen des Hauptsatzes 
bzw. des Zusatzes giiltig. Der Versuchsvektor u = (4, --., u,) habe genau k — 1 
bzw. n — k Zeichenwechsel. Aus u und dem transformierten Vektor 


At ==" (Wye a) 


1) G. Pétya und G. Szec6, Aufgaben und Lehrsdtze aus der Analysis, Bd. II (Springer, Berlin 
1925 [Grundlehren der mathematischen Wissenschaften]), siehe Abschnitt V, Kapitel 1, Aufgaben 7 


und 15. (Siehe auch Aufgabe 13.) 
2) L. Corrarz, Eigenwertaujgaben mit technischen Anwendungen ( 


1949), S. 289. 


Geest und Portig KG., Leipzig 
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bilde man die Quotienten 
Pays Li da 
a= u, da 


oat 


Sind saimtliche q, positiv, so liegt die k-te charakteristische Zahl zwischen Max(q;) 
und Min(q,). 

Auch den Collatzschen Einschliessungssatz fiir Differentialgleichungen!) kann 
man mit Hilfe der klassischen analytischen Oszillationssatze in analoger Weise 
modifizieren. Bei den anderen bekanntesten Einschliessungsverfahren, nament- 
lich bei denen, die sich auf den Rayleighschen Quotienten beziehen, ist jedoch 
eine solche Modifikation nicht ohne weiteres méglich. 


Summary 


Conditions on a linear algebraic eigenvalue problem are given, under which 
there are exactly k —1 changes of sign in the sequence of components of the 
k-th eigenvector. This is analogous to the oscillation theorems of differential 
equations. A class of difference equations which satisfies these conditions is 
defined. Finally a modification of a method of CoLLaTz is given, by means of 
which upper and lower bounds for the k-th eigenvalue may be derived from a 
trial vector having k — 1 sign changes in the sequence of its components. This 
paper is merely a summary of results; no proofs are given. 


(Eingegangen: 29. September 1953.) 
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Jahrestagung der Gesellschaft fiir angewandte Mathematik und Mechanik 
vom 20. bis 24. April 1954, in Miinchen 


Die Gesellschaft fiir angewandte Mathematik und Mechanik (GaMM) wird 
ihre nachste Jahrestagung in Miinchen in der Zeit vom 20. bis 24. April 1954 
veranstalten, Die 6rtliche Tagungsleitung wird in Handen von Prof. Dr. J. HEIn- 
HOLD (Technische Hochschule Miinchen) liegen. H. G6RTLER 


Internationale Tagung tiber Halbleiter in Amsterdam 
vom 29. Juni bis 3. Juli 1954 


Die Niederlandische Physikalische Gesellschaft, mit Unterstiitzung der Union 
internationale de Physique pure et appliquée und der Unesco, hat die Absicht, 
vom 29. Juni bis 3. Juli 1954 in Amsterdam eine Internationale Tagung iiber 
Halbleitey durchzufiihren. 

Die folgenden Forscher (von denen zwei unter Vorbehalt) haben zugesagt, 
Hauptvortrage zu iibernehmen: J. BARDEEN, W.H. Bratrain, H. B. G. Ca- 
sImIR, F. A. KROGER, D. PoLpER, M. ScH6n, W. SHOCKLEY, R. A. SMITH, H. J- 
Vink, H. WELKER. Hauptthemen der Tagung sind: Volumenrekombination; 
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Oberflachenleitfahigkeit; Oberflachenhaftstellen; Oberflachenrekombination : 
intermetallische Verbindungen; das Bandermodell fiir polare und nichtpolare 
Halbleiter; Photoleitfahigkeit in Halbleitern wie PbS; Pble, PbSe, ZnS, Cds; 
die Anwendung allgemeiner physikalischer und chemischer Gesetze auf die Her- 
_stellung von Halbleitern mit spezifischen Eigenschaften. 
; Ferner sind Fachausschuss-Sitzungen vorgesehen zur naheren Erérterung von 
_Einzelheiten. In diesen Sitzungen k6nnen auch kleinere Referate (15 Minuten) 
aus den verschiedenen Gebieten gehalten werden. Im weitern sind einige Ausfliige 
'beabsichtigt. 

Interessenten, welche dieser Tagung beiwohnen méchten oder einen Beitrag 
zu geben wiinschen, werden gebeten, sich so bald wie méglich an den Tagungs- 
sekretar (Dr. H. J. Vinx, Floralaan 142, Eindhoven, Holland) zu wenden. 

: H. J. Vink 
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Errata 


Experimentelle Untersuchung von Vereisungserscheinungen. Von 
DOMENIC MELCHER (ZAMP 2, Fasc. 6, 421 [1951)). 


Auf Seite 432 wird die wahrend einer Vereisung an einem Kabel auftretende 
- TemperaturerhGhung berechnet. 
Der Warmeentzug durch den Luftstrom ergibt sich rechnerisch zu 


Om af Aden G9) Oca call cus: 


wobei fiir AT die ohne Warmeabfluss pro Sekunde zu erwartende Temperatur- 
erhéhung eingesetzt wurde. Bei der Umrechnung der Warmeiibergangszahl ist 
leider ein Fehler im Betrag von zwei Zehnerpotenzen unterlaufen. 
: Q, fallt also hundertmal kleiner aus, als angegeben, und wird nicht von glei- 
cher Gréssenordnung wie die frei werdende Kristallisationswarme Q,x. Damit 
wird im stationaéren Fall die gesuchte Temperaturerhéhung richtig von der 
_ Gréssenordnung einiger Zehntelgrad. 
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Statik der rahmenartigen Tragwerke. Von J. PIRLET (Springer-Verlag, 
Berlin 1951). 168 S., 80 Abb.; DM 24. 
Der Verfasser behandelt die Berechnung von «rahmenartigen Tragwerken», 
die sich aus elastisch eingespannten Balken oder Staében zusammensetzen und 
deren Knotenpunkte wohl elastisch verdrehbar, aber unverschieblich sind. Die 
-Berechnung dieser hochgradig statisch unbestimmten Systeme wird schrittweise 
durchgefiihrt; wird nur ein einziges Balkenfeld als belastet betrachtet, so ergeben 
sich seine beiden Einspannmomente aus den « Einspannungsgraden» der beiden 
Stabenden, die ihrerseits schrittweise, von den Auflagerpunkten ausgehend, be- 
-rechnet werden kénnen. In umgekehrter Richtung kann anschliessend die « Fort- 
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pflanzung» der Momente gefunden werden. Aus dem Umstand, dass die Ein- 
spanngrade eines Feldes durch kleine Anderungen der Einspanngrade der an- 
schliessenden Stabe nur unwesentlich beeinflusst werden, lasst sich ein verein- 
fachtes Berechnungsverfahren durch angenaherte Bestimmung der Einspann- 
grade aufstellen. Das Erfreuliche an diesem Buche liegt darin, dass es auf den 
Grundlagen der klassischen Baustatik aufbauend zu einer klaren, tibersichtlichen 
und strengen Berechnung der untersuchten besonderen Tragwerksform Per 

. Stusst 


Elasticity in Engineering. By Ernest E. SECHLER (John Wiley & Sons, 
New York, 1952). 419 Seiten, mit vielen Abbildungen; $8.50. 


Dieses Buch — es ist das neueste aus der CALCIT-Serie (Guggenheim Aeronau- 
tical Laboratory, California Institute of Technology) — ist fiir Ingenieure be- 
stimmt, die auf allen Gebieten der héheren Festigkeitslehre arbeiten. 

Das Buch besteht aus drei Teilen. Im ersten Teil werden die gebrauchlichen 
Beziehungen in verschiedenen Koordinatensystemen hergeleitet, welche fiir zwei- 
und dreidimensionale Probleme das Gleichgewicht, die Vertraglichkeit und die 
elastische Energie ausdriicken. Ferner wird der Begriff der Spannungsfunktion 
eingefiihrt. Der zweite Teil befasst sich mit Systemen, die sich im stabilen Gleich- 
gewicht befinden. Hier werden zweidimensionale Probleme, Balkenbiegung, 
statisch unbestimmte Fachwerke, Torsion von Wellen, behandelt und die Grund- 
gleichungen der Platten- und Schalentheorie diskutiert. Der dritte Teil umfasst 
Instabilitatsprobleme bei gedriickten Saulen, gedriickten und gebogenen Balken 
sowie bei Platten und Schalen. 

Das Buch setzt offensichtlich wenig Vorbildung voraus. Obschon die Auswahl 
des Stoffes im allgemeinen dem Zweck entspricht, scheint es dem Referenten, 
dass gewisse elementare Dinge, wie der Mohrsche Satz iiber die elastische Linie 
und die Lésung einfacher statisch unbestimmter Probleme, eher in ein Buch iiber 
‘elementare Festigkeitslehre gehéren. 

Das Buch ist dadurch gekennzeichnet, dass die Zwischenrechnungen bei der 
Aufstellung einer Gleichung oder bei der Lésung eines Beispiels stets in extenso 
aufgefiihrt werden. Andererseits ware es zu begriissen, wenn der begleitende Text 
die physikalischen Aspekte besser herausarbeiten wiirde. Die am Ende jedes 
Kapitels zusammengestellten Ubungsaufgaben sind dagegen dusserst wertvoll. 

E. P. Popov 


Priifungs- und Ubungsaufgaben aus der Mechanik des Punktes und 
des starren Kérpers. Von K. FEDERHOFER (Springer-Verlag, Wien 1950/51). 
I. Teil: Statik, 130 S., 243 Abb.; sFr. 10.—. II. Teil: Kinematik und Kinetik 
des Punktes, 103 S., 105 Abb.; sFr. 10.-. III. Teil: Kinematik und Kinetik starver 
Systeme, 139 S., 191 Abb.; sFr. 10.-. 


, 


Die vorliegende Aufgabensammlung zur Mechanik ist in drei getrennt erschie- 
nene Teile gegliedert, von denen der erste der Statik, der zweite der Kinematik 
und Kinetik des Punktes und der dritte der Kinematik und Kinetik starrer 
Systeme gewidmet ist. Jeder Teil enthalt getrennt die Aufgaben und ihre Lésun- 
gen, wobei den letzteren eingehende Erlauterungen iiber den zweckmassigen 
Loésungsweg, die stellenweise den Charakter eigentlicher Darstellungen der be- 
handelten Probleme annehmen, vorangestellt werden. Die erfolgreiche Beherr- 
schung und Anwendung der Lehren der Mechanik setzt eine selbstandige Erar- 
beitung des Stoffes voraus; dazu kann die Beispielsammlung K. FEDERHOFERS, 


j 
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_ die im Laufe von drei Jahrzehnten erfolgreicher Lehrtatigkeit entstanden ist, 
_ eine wertvolle Hilfe sein. F.. Stiissi 


Grundlagen der Elektroakustik. Von W. REICHARDT (Akademische Ver- 
lagsgesellschaft Geest & Portig KG., Leipzig 1952). 464 S., 311 Abb.; DM 32.-. 


Das vorliegende umfangreiche Werk ist ein eigentliches Handbuch der 
Elektroakustik und bietet mehr als nur eine Einfiihrung in die Grundlagen dieses 
Gebietes, wie im Vorwort allzu bescheiden bemerkt ist. Dies zeigt schon das um- 
fangreiche Literaturverzeichnis mit iiber 400 Hinweisen, das allerdings keinen 
Anspruch auf Vollstandigkeit machen kann und will. 

Der erste Teil des Buches behandelt die Grundlagen der physikalischen Aku- 
stik. Gleich zu Anfang wird vorgeschlagen, den allgemein bekannten und ge- 
brauchlichen Begriff der akustischen Impedanz durch den Ausdruck «Schall- 
wellenstandwert» zu ersetzen; die selten verwendete akustische Admittanz wird 
«Schallwellenmitgang» genannt. Die Wiinschbarkeit, eine international einge- 
fuhrte und bewahrte Terminologie durch rein deutsche neue Ausdriicke zu er- 
- setzen, kann fiiglich bezweifelt werden. Ferner ist der Versuch, die Schalldam- 
mung zusammen mit Reflexion, Beugung, Brechung und Ausbreitung in der 
4 Atmosphare zu behandeln, nicht nur aus wissenschaftlichen, sondern vor allem 
_ auch aus didaktischen Gesichtspunkten abzulehnen. Der Zusammenhang des 
_ Kapitels « Nichtlinearitaten», das den ersten Teil des Buches abschliesst, mit der 
_ ubrigen Materie ist schwer einzusehen. 

Als weiterer Hauptabschnitt folgt «Das Schallempfinden des Menschen», in 
dem sich schone Darstellungen tiber den Aufbau und die wichtigsten Eigenschaf- 
_ ten des Ohres finden. Trotz dem recht beachtlichen Umfang des Kapitels «Laut- 
_ starke» vermisst man aber Angaben iiber die praktisch wichtige Lautstarke- 
_ berechnung sowie besonders auch iiber Verstandlichkeitsberechnungen, die heute 

in der Technik mehr und mehr Eingang finden. Die Darstellung der resultieren- 
_ den Lautstarke mehrerer Schallquellen ist leider bei der viel zu primitiven Vor- 
_ stellung einer energiemassigen Addition stehengeblieben. 
Der vierte Teil des Buches ist den Schallwandlern gewidmet, wobei zuerst die 
_ Kolbenmembran, insbesondere ihre Richtwirkung zur Sprache kommt. Ein um- 
fangreiches Kapitel, «Mechanische Schwingungsgebilde», behandelt im wesent- 
_ lichen nur den schwingenden Massenpunkt, wobei die mechanisch-elektrischen 
_ Analogien sehr breit diskutiert werden. Leider fehlt die Behandlung von kon- 
 tinuierlichen Schwingungsgebilden, wie zum Beispiel die gespannte Membran 
und die eingeklemmte Platte, die im Wandlerbau wichtige Elemente sind. An- 
_ schliessend folgen die Kontaktmikrophone, das Telephon, der elektromagnetische 
_ Lautsprecher, die elektrodynamischen Mikrophone in ihren verschiedenen be- 
_ kannten Varianten, der elektrodynamische Lautsprecher, das Kondensatormikro- 
_ phon und schliesslich das Kristallmikrophon, eine Systematik der Behandlung, 
_ die nicht immer gliicklich ist. 

Im Abschnitt «Mikrophoneinsatz» vermittelt der Verfasser offensichtlich 
_eigene wertvolle Erfahrungen; ganz am Rande wird hier auch die statistische 
- Raumakustik noch gestreift. Auch das folgende Kapitel, « Beschallungstechnik», 
_ bringt viel wertvolles Material, wobei besonders die in den letzten Jahren immer 
_haufiger verwendeten Strahlergruppen erwahnt seien. Schade, dass hier nicht 
auch der sogenannte «Haas-Effekt» aufgefiihrt wird, der eine der wichtigsten 
-neueren Erfahnungen auf diesem Gebiete ist. 

Das Buch schliesst mit der Behandlung der Schallaufzeichnung, wobei das 
- Schwergewicht richtigerweise auf dem Grammophon und dem Magnetophon liegt. 
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Im ganzen Buche wird konsequent das Giorgi-MaBsystem angewendet, tiber 
das eine vorziigliche klare und kurze Einleitung vorhanden ist. Der schon ein- 
gangs erwahnte Dualismus zwischen Handbuch einerseits und elementaren Grund- 
lagen andererseits muss zu einigen Widerspriichen fiihren, die die Lekttre des ; 
Buches erschweren. Seine didaktischen Mangel lassen es als Lehr- und Studien- 
mittel weniger geeignet erscheinen, dagegen wird es dem in der Praxis stehenden 
Ingenieur und Techniker ein sehr willkommenes Hilfsmittel sein und ihm wert- — 
volle Unterlagen und Anregungen vermitteln. W. Furrer 


Dynamik des Bogentragers und Kreisringes. Von K. FEDERHOFER 
(Springer-Verlag, Wien 1950). 179 S., 35 Abb.; sFr. 23.50. 


Der Verfasser unternimmt es hier, zu untersuchen, wie sich die Besonderheiten 
der Biegung und Verdrehung, wie sie bei unsymmetrischen offenen Querschnitten 
auftreten und die durch den Begriff des Schubmittelpunktes charakterisiert wer- 
den kénnen, auf die Eigenschwingungen gekriimmter Stabe auswirken. Es liegt 
in den Schwierigkeiten der Integration der bei der mathematischen Formulierung 
des Problems auftretenden Differentialgleichungen begriindet, dass die Unter- 
suchung von einschrankenden Voraussetzungen (kreisf6rmige Bogenachse, kon- 
stanter Querschnitt) ausgehen muss; die Einfliisse von Abweichungen gegentiber 
diesen Voraussetzungen kénnen erganzend untersucht werden. Das Buch FEDER- 
HOFERS stellt eine willkommene erste Verbindung zwischen der verallgemeinerten 
Biegungslehre und der Dynamik des gekriimmten Stabes her. F, Sttisst 


Theoretical Nuclear Physics. By J. M. BLatr and V. F. WEIssKoprF (John 
Wiley & Sons, New York; Chapman & Hall, London, 1952). 864 pp.; $12.50. 


Von einem grundsatzlichen Verstandnis der im Kern wirkenden Krafte sind 
wir noch weit entfernt, und eine Deutung der Phanomene ist deshalb nicht, wie 
etwa in der Atomphysik, ein rein mathematisches Problem auf Grund wohldefi- | 
nierter physikalischer Vorstellungen. Vielmehr miissen fiir jeden Problemkreis 
eigene Modelle erdacht werden, die eine méglichst weitgehende Beschreibung 
der betreffenden Phanomene erlauben, deren Zusammenhang jedoch oft dunkel 
ist und einer spateren Synthese iiberlassen bleiben muss. Aus diesen die bewahrten 
und hoffnungsvollen herauszuschalen und in eine einheitliche Darstellung zu 
giessen, ist ein Unterfangen, das nur bei tiefgehender Kenntnis des Gebietes er- 
folgreich sein kann. Den Verfassern des vorliegenden Buches ist das véllig ge- 
lungen; sie haben ein Werk geschaffen, das sowohl als Lehrbuch fiir Studierende 
wie als Nachschlagewerk fiir experimentelle und theoretische Physiker einzig 
dasteht. 

Das Buch gliedert sich im wesentlichen in vier Teile. In einem ersten Teil 
(Kapitel I bis VII und XIV) wird, in induktivem Aufbau von einfacheren zu 
komplizierteren Systemen, versucht, méglichst viel Information iiber die Kern- 
krafte zu gewinnen. Theoretische Spekulationen iiber die Natur dieser Krafte 
(im Sinne der Mesontheorie) fehlen, da der heutige Stand der Feldtheorien es 
nicht gestattet, irgendwelche quantitativen Folgerungen aus solchen zu ziehen. 
Konsequenterweise ist auch eine Diskussion der experimentellen Tatsachen der 
Mesonenphysik unterblieben; das Buch beschrankt sich auf das Gebiet kleiner 
Energien (< 50 MeV), mit Ausnahme eines Kapitels iiber die Nukleonstreuung 
bei hohen Energien, welche eines der Hauptratsel der Kernphysik bildet. Ein 
zweiter Teil (Kapitel VIII bis XI) behandelt die allgemeine Theorie der Kern- 
reaktionen. Ausser einigen exakten Theoremen handelt es sich hierbei vor allem 
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um eine phanomenologische Beschreibung der Tatsachen auf Grund einiger fun- 
damentaler Annahmen (Compoundkern). Ein dritter Teil (Kapitel XII) enthalt 
die Theorie der Wechselwirkung von Kernen mit dem elektromagnetischen Strah- 
lungsfeld (y-Strahlen). Dieser Abschnitt ist besonders bemerkenswert, da er viel 


_ bisher unpubliziertes Material enthalt, welches zum Teil auch Irrtiimer 4lterer 


Arbeiten klarstellt. Ausserdem zeichnet er sich durch einen besonders eleganten 
Formalismus zur Beschreibung von Multipolstrahlung aus. Der vierte Teil end- 
lich (Kapitel XIII) behandelt den f-Zerfall. Auch hier ist viel Konfusion in der 
alteren Literatur vermieden durch konsequente Unterscheidung der «verbotenen » 
und der «relativistischen» Ubergange. 

Das ganze Werk ist leicht zu lesen, und trotz weitgehend strenger mathema- 
tischer Behandlung werden dem Leser lange formale Rechnungen erspart, 
eventuell durch Vorfiihren vereinfachter, aber typischer Beispiele. (Ob allerdings 
die weitere Tendenz, jedes formale Resultat durch eine «anschauliche» Plausibili- 
tatsbetrachtung nochmals herzuleiten, nicht etwas tibertrieben wird, mag Ge- 
schmacksache sein). Ausfiihrliche Tabellen der verwendeten Symbole erleichtern 


_ das Nachschlagen sehr. Das Buch wird wohl auf Jahre das Standardwerk tiber 


theoretische Kernphysik bleiben. R. Schafroth 


Physik diinner Schichten, Teil I: Herstellung, Dickenmessung und optische 
Eigenschaften. Von HERBERT Mayer (Wissenschaftliche Verlagsgesellschaft 
MbH., Stuttgart 1950). 326 S., 142 Abb.; sFr. 37.50. j 

Diinne Schichten spielen heute auf vielen Gebieten der Technik eine bedeu- 
tende Rolle. Als Beispiel seien die reflexvermindernden Schichten in der Optik 
genannt. Die Vergiitung von Objektiven ist in den letzten Jahren zur Selbstver- 
standlichkeit geworden. Insbesondere bei Objektiven mit grésserer Linsenzahl 
ist diese Vergiitung von grosser Bedeutung, indem sie einen betrachtlichen Licht- 
gewinn mit sich bringt. 

Im vorliegenden Werk, Band 4 der Sammlung « Physik und Technik», wird 
erstmalig eine zusammenfassende Darstellung der Physik der diinnen Schichten 
gegeben. 

Nach einem Uberblick iiber die historische Entwicklung und einer Bestim- 
mung des Begriffes «diinne Schicht» wird im zweiten Kapitel eine Ubersicht iiber 
die wichtigsten Verfahren zur Herstellung dieser Schichten gegeben: Aufdampfen, 
Kathodenzerstaubung, Adsorption usw. 

Im nachsten Kapitel werden die Verfahren zur Messung der Dicke dinner 
Schichten behandelt: mechanische Messung, elektrische Methoden, chemische 
Verfahren und optische Methoden. Die letzten Methoden umfassen dabei: Inter- 
ferenzfarben-, Interferenzstreifen-, Intensitats- und Polarisationsmethode. Im 
gleichen Kapitel wird ferner die wichtige Grosse der Oberflachenrauhigkeit be- 
handelt. 

Im vierten Kapitel werden dann die Grundlagen der Optik diinner Schichten 
dargestellt und mit den experimentellen Ergebnissen verglichen: homogene 
Schichten, ihre Reflexions- und Absorptionsvermogen sowie ihre Durchlassigkeit, 
inhomogene und anisotrope Schichten. Es werden die Methoden zur Bestimmung 
der optischen Konstanten dieser Schichten dargelegt. Schliesslich geht der Autor 
auf Anwendung diinner Schichten in der Optik ein: Erhéhung und Verminde- 
rung der Reflexion durch diinne Schichten, diinne Schichten als Interferenz- 
filter sowie als Polarisatoren. 

Das Werk ist mit 142 Abbildungen reich illustriert. Es besitzt einen Anhang 
mit 20 wichtigen Tabellen tiber Reflexionsvermégen, Absorptions- und Brechungs- 
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index von hdufig verwendeten Metallen. Das Literaturverzeichnis mit tiber 600 
Hinweisen auf Arbeiten auf dem Gebiete diinner Schichten ist sehr vollstandig. — 

Im vorliegenden Werk ist eine klare Darstellung iiber das wichtige Gebiet _ 
der diinnen Schichten gegeben. Wir méchten es jedem, der auf diesem Gebiete 
tatig ist oder sich in dieses einarbeiten will, voll empfehlen. N. Schaetts | 


Zur Struktur und Materie der Festkérper. Diskussionstagung der Sek- 
tion fiir Kristallkunde der Deutschen Mineralogischen Gesellschaft (Springer-_ 
Verlag, Berlin 1952). 304 S., 95 Abb.; DM 28.60. | 


Die an der Diskussionstagung der Sektion fiir Metallkunde der Deutschen 
Mineralogischen Gesellschaft am 1. und 2. Mai 1951 in Frankfurt am Main ge- | 
haltenen Vortrage sowie die zugehérigen Diskussionsbeitrage sind in Buchform 
herausgekommen. Ahnlich wie der Band Der feste Kérper der Ziircher Tagung — 
(1938) zeigt dieses Buch einen Querschnitt durch den heutigen Stand der Kennt- 
nisse und Probleme der reinen und angewandten Kristallstrukturforschung. Bei- 
trage zur allgemeinen Symmetrie- und Strukturlehre lieferten P. Nicci (Ein- 
fiihrung von Charakterentafeln fiir Fourier- und Patterson-Analyse), C. HERMANN 
(Translationsgruppen in n Dimensionen) und R. HosEMANN (Versuch einer allge- 
meinen Interferenztheorie, speziell im Hinblick auf Parakristalle, hochmolekulare 
und kolloide Stoffe). Die Wachstumstheorie wurde durch W. KossEL beleuchtet 
(Beobachtung von Kristallkugeln als Forschungsmittel fiir Oberflachenvorgange) und 
durch I. N. Stransk1 und Mitarbeiter aufs wertvollste erganzt mit einer Arbeit 
iiber die Verdampfung von Kristallen. Ins}Gebiet der Realbauforschung fallt der 
Beitrag von H. JaGopzinsx1, der die kooperative Fehlordnung in Kristallen 
theoretisch beleuchtete. Aus der physikalisch-chemischen Kristallforschung be- 
richtete einerseits H. WittE (Einfluss der Brillouin-Zonen auf die Eigenschaften 
von Legierungen), in der Diskussion durch K. SCHUBERT noch wesentlich erganzt, 
und anderseits A. SMEKAL, der die aus der Kristallchemie stammenden Erkennt- 
nisse liber die Bindungsverhaltnisse auf die Glaser anwandte und versuchte, eine 
bindungsmassige Theorie des Glaszustand-Vorkommens zu geben. Diese kurzen 
Hinweise auf den Inhalt der Einzelkapitel mégen geniigen, um zu zeigen, dass hier 
eine wertvolle Zusammenstellung modernster kristallographisch-physikalisch-che- 
mischer Arbeiten vorliegt, die jedem mit Kristallen arbeitenden Chemiker, Phy- 
siker und Mineralogen gewinnbringend sein werden. W. Epprecht 


Introductory Nuclear Physics. By D. Hatiipay (J. Wiley & Sons, Inc., 
New York, 1950). 558 pp.; $6.50. 


Dieses Buch stellt eine klare und leichtverstaéndliche Einfiihrung in die 
moderne Kernphysik dar. Es enthalt in knapper und einfacher Formulierung 
sowohl die wichtigsten theoretischen Vorstellungen, die man sich heute iiber den 
Aufbau der Kerne sowie den Mechanismus der Kernumwandlungen macht, als 
auch eine Beschreibung der grundlegenden Experimente, welche zu diesen Vor- 
stellungen gefiihrt haben. Auch die fiir ein Verstandnis der Kernvorgange uner- 
lasslichen Nachbargebiete, wie Wechselwirkung zwischen Strahlung und Materie, 
sowie kosmische Strahlen kommen kurz zur Behandlung. Zwei weitere Abschnitte 
befassen sich mit den experimentellen Methoden zum Nachweis geladener Partikel 
und Quanten sowie der kiinstlichen Erzeugung energiereicher Teilchen. Der Leser — 
als solcher kommt fiir dieses Buch der Studierende mittlerer Semester wie auch 
der in der Praxis stehende Ingenieur mit physikalischen Interessen in Frage — er- 
fahrt damit alles Wesentliche, das er fiir eine erste Orientierung wie auch fiir ein 
nachfolgendes detailliertes Studium auf diesem Gebiet benotigt. H. Waffler 
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A. Magnetisierungserscheinungen — B. Streuungserscheinungen 
C. Uberspannungserscheinungen — D. Betriebseigenschaften 
E. Sonderschaltungen — F. Spannungsregelung — G. Experi- 
mentelle Untersuchung — H. Aufbau des Transformators 
I. Kiithlung und Erwarmung — K. Beanspruchung und Wahl 


der Isolierung — L. Entwurf. 
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VERLAG BIRKHAUSER - BASEL/STUTTGART 


Wissenschaftliche Zeitschriften 
aus dem Verlag Birkhauser, Basel/ Stuttgart 


Technica, Illustrierte technische Zeitschrift. Erscheint vierzehntiglich. Abonne- 
mentspreis pto Jahr Fr. 14.40, Ausland Fr. 18.—. Einzelnummer 60 Rp., 
Ausland 75 Rp. . 


Helvetica Physica Acta. Edita a Societate Physica Helvetica. Erscheinen seit ~ 
1928. 23 X15,5 cm. Abonnementspreis pro Jahr Fr. 44.-; Ausland 
Fr. 52.-; Einzelnummer Fr. 11.— 


Archiv der Mathematik. Herausgegeben von A. Kneser und W. Suis. Exscheint 
sechsmal jahrlich, seit 1948. 24 x 17 cm. Abonnementspreis noes Jahr 
Fr, 60.-; Binzelnummer Fr. 12.— 


Elemente der Mathematik. Revue de mathématiques élémentaires. Rivista di 
matematica elementare. 24 x 16,5 cm. Erscheinen alle 2 Monate. Jahr- 
lich Fr. 10.-, Ausland Fr. 14.-; Einzelnummer Fr. 2.50. 


Experientia. Revue mensuelle des sciences putes et appliquées. Monatsschrift fur 
das gesamte Gebiet der Naturwissenschaft. Rivista mensile di scienze 
pure e applicate. Monthly Journal of Pure and Applied Science. Edi- 
totes: A. v. Muralt, Bern — L. Ruzicka, Zurich — J. Weigle, Genéve. 
Erscheint seit 1945. Abonnementspreis pro Jahr Fr. 32.-, Ausland 
Fr, 38.—, Einzelnummer Fr. 3.50. 


Schweizerische paléontologische Abhandlungen — Mémoires suisses de 
paléontologie. Herausgegeben von der Schweizerischen Paldontolo- 
gischen Gesellschaft. Erscheinen seit dem Jahre 1874, pro Jahr durch- 
schnittlich ein Band. 32 x 24 cm. Preis pro Band Fr. 75.—. 


Eclogae Geologicae Helvetiae. Herausgegeben von der Schweizerischen Geo- 
logischen Gesellschaft. Erscheinen seit 1888. 24,5 x17 cm. Jahrlich 
: 1 Vol. mit 2 Heften 4 ca. Fr. 20.— (je nach Umfang). 


Schweizerische Zeitschrift fiir Hydrologie. Revue suisse d’Hydrologie. Her- 
ausgegeben von der Hydrobiologischen Kommission der Schwei- 
zerischen Natutforschenden Gesellschaft. Jahrlich erscheinen 2 Hefte, 
seit 1920. Abonnementspreis: Inland Fr, 24.-, Ausland Fr. 28.-. 
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Werke der reinen nd angewandten Mathematik und Physik 


Mathematik 


BieBrRBACH, L.: Theorie der geometrischen Konstruktionen . 
BLASCHKE, W.: Projektive Geometrie ..--.+-++-++++ ++ 5 


Locuer-Ernst,L.: Differential- und Integralrechnung . ... . 
Locuer-Ernst,L.: Einfiihrung in die freie Geometrie ebener Kurven 
Ostrowsk1,A.: Vorlesungen tiber Diff.- u. Integralrechng., Band I 
(Band III erscheint demnachst) Band II 

PUNE L.: Gesammelte mathematische Abhandlungen, Band I 
(Band Ill in Vorbereitung) BandII 


_STIEFEL, E.: Lehrbuch der darstellenden Geometrie. .. . . 
‘Trost, E.: Primzahlen ...... ree Sia abe Newiawcs Se ences Fr. 13.50 


VoELKER, D., Dorrscu, G.: Die zweidimensionale legate Eek 
Zwinaol,E.: Versicherungsmathematik ...........-. 


Physik und Astronomie 


GomsAs,P.: Theorie und Lésungsmethoden des Mehrteilchenpro- 
blems der Wellenmechantk <7. .) seu. <1 ssn toga Cen 
NIETHAMMER, TH.: Astronomisch-geogr. Ortsbestimmung . . . - 


geb. Fr. 18.70 
geb. Fr. 19.60 


Buascuxke, W.: Analytische Geometrie. ..--.++++-+ = geb. Fr. 19.60 
CaratHtopory, C,: Funktionentheorie, Bd. I geb. Fr. 39.50, Bd. II geb. Fr. 27.05 
‘Caratutopory, C.: Mass und Integral und ihre Algebrdisierung . erscheint 1954 
Dortscu,G.: Handbuch der Laplace-Transformation. Band I: 
Die theoretischen Grundlagen der Laplace-Transformation . . geb. Fr. 83.20 
Finser, P.: Kurven und Flachen in allgemeinen Raumen ... geb. Fr. 15.40 
' Furrer,R.: Analytische Geometrie der Ebene und des Raumes . geb. Fr. 24.95 
Herrmann, H.: Ubungen zur projektiven Geometrie ..... geb. Fr, 27— 
Hurwitz,A.: Mathematische Werke,2 Bande ....... - eb: je Fr. 54.10 
LinveER,A.: Statistische Methoden, 2. Auflage 1951. ...... geb. Fr. 31.20 


geb. Fr. 49,90 
geb. Fr. 12.50 
geb. Fr. 36.40 
geb. Fr. 69.70 
geb. Fr. 56.15 
geb. Fr. 56.15 
geb. Fr. 29.65 


geb. Fr. 47.85 
geb. Fr. 31.20 


geb. Fr. 33.30 
geb. Fr. 36.40 


VON DER PAHLEN, E.; Einfithrung in die Dynamik von Sternsystemen geb. Fr. 39.50 


Wavpmeirr, M.; Einfihrung in die Astrophysik ........- 
Kongressbericht der Internat. Konferenz iiber Kernphysik und 
theoretische Physik, Basel 1949. Broschiert 


~ Ingenteurwissenschaften 


Bopea,E.: Giorgis rationales MKS-MaBsystem. Broschiert. . . . 
GOLDSTEIN, J.: Die Messwandler . ........ oo ip Srey ate te 
JAEGER, Cu.: Technische Hydraulik. .. 5.) 3. ee 
KuUHLMANN, K.: Theoretische Elektrotechnik. Band II: Grundziige 


der Theorie der Wechselstréme u. des einphasigen Transformators 


Band III: Grundziige der Theorie elektrischer Maschinen. . . 
MEISSNER, E., und Z1EGLER, H.: Mechanik. Band I: Statik der star- 
ren, fliissigen und elastischen Kérper. . ..........- 
Band II: Dynamik der starren Kérper. .......... 
Band III: Dynamik der Systeme . .......... mete 
MicuazL, W.: Ortskurvengeometrie in der komplexen Zahlenebene 
SANGER,R,: Ballistische Storungstheorie ........... 
Stuss1,R.: Vorlesungen iiber Batstatik, BandI ....... 
RUSTERHOLZ, A,: Elektronenoptik, Band I: Grundziige ay yee eee 
RicnTer, R.: Elektrische Maschinen Band I 
Band II geb. Fr. 46.80, Band III 
Vipmar, M.; Die Gestalt der elektrischen Freileitung ..... . 
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geb. Fr. 49.40 


Fr. 16.65 


Fr. 28.10 
geb. Fr. 29.10 
geb. Fr. 52.— 


geb. Fr. 27.55 
geb. Fr. 74.90 


geb. Fr. 37.45 
geb. Fr. 33.80 
geb. Fr. 46:80 
geb. Fr. 13.50 
geb. Fr. 15.10 
geb. Fr. 40.05 
geb. Fr. 31.20: 
geb. Fr. 49.40 
geb. Fr. 26.— 
geb. Fr. 19.75 
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